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FOREWORD 

Tills  technical  report  is  being  issued  under  the  Control  Systems 
Development  Branch,  Flight  Control  Division,  Air  Force  Flight  Dynamics 
Laboratory,  Air  Force  Wright  Aeronautical  Laboratories,  Wright-Patterson 
Air  Force  Base,  Ohio,  under  Project  2A03,  "Flight  Control,"  Task  2A0302, 
"Flight  Control  Systems  Development."  This  report  was  a dissertation 
prepared  by  Ernest  F.  Moore,  M.  S.,  Lt  Colonel,  USAF,  and  presented  to  the 
Faculty  of  the  School  of  Engineering  of  the  Air  Force  Institute  of  Technology 
in  partial  fulfillment  of  the  requirements  for  the  degree  of  Doctor  of 
Philosophy. 

This  dissertation  is  part  of  a sequence  of  studies  being  conducted 
at  the  Air  Force  Institute  of  Technology  to  build  a sound  technology 
base  concerning  the  propagation  of  signals  in  fluid  transmission  lines. 

The  central  objective  of  this  dissertation  was  to  develop  the  frequency 
and  transient  response  of  transmission  lines  of  arbitrary  cross  section 
carrying  fully  developed  laminar  or  turbulent  mean  flow  of  a compressible 
or  incompressible  fluid. 

A legacy  of  many  computer  programs  and  large  quantities  of  data, 
much  too  voluminous  for  inclusion  in  this  work,  has  been  organized  and 
bound  for  the  convenience  and  use  of  the  next  AFIT  student  who  desires 
to  continue  tills  endeavor.  It  is  safe  to  say  that  the  solution  of  one 
problem  surfaces  at  least  two  more,  so  that  the  quest  is  not  ended  with 
this  work. 

There  are  many  people  to  thank.  Captains  G.  R.  Farney  and  J.  D. 
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members  of  the  AFIT  Research  Advisory  Committee,  composed  of  Dr.  M.  E. 
Franke,  Dr.  A.  J.  Shine,  Dr.  C.  H.  Houpis,  Major  J.  T.  Karam,  Jr.,  and 
lLt.  R.  Merz  for  their  guidance.  Special  thanks  are  extended  to  the 
chairman,  Dr.  Franke,  for  his  patience  and  encouragement,  and  for  the 
sound  advice  given  throughout  this  work. 
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laminar  flow  condition 
low  frequency  or  long  time 
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annular  line 

circular  or  outside  radius 
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receiving  station  in  line 
sending  station  in  line 
lossless  condition 
turbulent  flow  condition 
thermal  or  temperature 
arbitrary  line 
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( ) average  or  normalized 

instantaneous  deviation  from  average  or  alternate  form 
( )*  reference  conditions 

( )°  degrees  temperature 


THE  SMALL -SIGNAL  RESPONSE  OF  FLUID  TRANSMISSION  LINES 


INCLUDING  DEVELOPED  MEAN  FLOW  EFFECTS 

I.  Introduction 

Background 

The  propagation  of  small  signals  through  ducts  or  lines  filled 
with  a fluid  has  been  of  interest  to  scientists  and  engineers  for  some 
time.  With  the  advent  of  fluidics  and  fluerics  in  the  late  1950's  and 
up  to  the  present  time,  many  investigations  have  been  completed,  both 
analytically  and  experimentally,  in  order  to  predict  and  measure  the 
parameters  governing  small  signal  propagation. 

The  modeling  of  fluid-filled  transmission  lines  has  progressed 
from  the  basic  lossless  line  model,  to  the  line  whose  losses  might  be 
called  linear  with  mean  velocity,  to  the  present  frequency-dependent 
friction  models.  Solutions  have  been  obtained  in  the  frequency  domain 
for  lines  of  circular  and  rectangular  cross  section,  as  these  cross 
sections  are  commonly  found  in  fluid  systems.  For  the  fluid-filled 
circular  line,  the  transient  response  to  small  amplitude  impulse  and 
step  inputs  has  been  obtained  in  various  forms. 

Recently,  investigations  have  been  carried  out  and  models 
developed  to  determine  and  to  describe  the  response  of  liquid-filled 
circular  lines  to  oscillatory  inputs,  where  the  inputs  are  superimposed 
on  a steady  developed  flow,  and  where  the  steady  developed  flow  may  be 
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either  laminar  or  fully  turbulent.  A limited  amount  of  experimental 
data  has  provided  some  degree  of  credence  to  these  models . 


Problem  Description 

During  the  literature  search  preceding  this  work,  it  appeared 
that  a majority  of  the  investigations  were  very  specialized  and  that 
an  effort  to  synthesize  the  results  obtained  by  others  was  in  order. 

Much  of  the  previous  work  relating  to  transient  response  in  the  time 
domain,  and  to  combined  oscillatory  and  steady  developed  flow,  was 
restricted  to  liquid-filled  circular  lines. 

A solution  for  the  small-signal  response  of  fluid  lines  of 

j 

concentric  annular  cross  section  appeared  to  be  a useful  first  step  in 
synthesizing  and  expanding  the  applicability  of  transmission  line 
modeling.  Since  the  annular  line  includes  the  circular  line  as  a 
special  case  (when  the  inside  radius  is  equal  to  zero)  and  approaches 
the  parallel  plate  configuration  in  its  other  extreme,  a solution  was 
first  obtained  for  this  line. 

In  studying  the  solutions  available  for  circular  and  rectangular 
lines  as  well  as  the  solution  for  annular  lines  given  herein,  it  appeared 
that  a method  of  expressing  the  frequency  dependent  response  of  lines  of 
arbitrary  cross  section  in  terms  of  only  the  circular  line  might  be 
developed.  Through  the  use  of  the  various  analogies  between  electrical 
transmission  line  theory  and  fluid  dynamics,  and  in  particular  by 
constructing  and  analyzing  the  equivalent  electrical  circuit  of  a 
gene-ralized  fluid  transmission  line,  such  a method  was  found. 

It  followed  that  if  the  frequency  dependence  of  a line  of 
arbitrary  cross  section  could  be  expressed  in  terms  of  an  equivalent 
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circular  line,  then:  (1)  The  small  signal  transient  (time  domain) 
response  of  lines  of  arbitrary  noncircular  cross  section  might  be  found, 
since  this  response  is  known  for  circular  lines;  (2)  The  characteristic 
impedance  of  lines  of  noncircular  cross  section  may  be  expressed  in 
terms  of  equivalent  circular  lines;  consequently,  impedance  matching  of 
lines  of  different  cross  section  may  be  possible;  (3)  The  computer 
modeling  of  all  lines  of  any  cross  section  in  a fluid  system  might  be 
reduced  to  modeling  the  equivalent  circular  lines;  (4)  A method  of 
determining  the  small  signal  response  of  fluid  lines  of  arbitrary  cross 
section,  carrying  fully  developed  mean  flow  might  be  found. 

All  of  these  possibilities  were  investigated  and  the  associated 
problems  were  solved,  subject  to  the  restrictions  specified  herein.  It 
is  noted  that  where  possible  the  simplest  model  achievable  is  used,  in 
accord  with  Oldenberger' s observation  that  only  those  models  which  are 
reasonably  uncomplicated  find  wide  use.  This  observation  is  particularly 
appropriate  to  fluid  systems,  which  generally  contain  a myriad  of 
components  and  lines,  some  of  which  are  usually  noncircular  in  cross 
section. 

Organization  and  Contents 

Since  this  work  covers  a rather  broad  range  of  transmission  line 
problems,  specific  background  information  and  references  will  be 
introduced  as  appropriate.  Further,  where  specific  examples  are  deemed 
necessary,  these  will  be  given  in  terms  of  the  more  common  rectangular 
and  annular  cross  sections  for  which  a relatively  larger  body  of 
technical  reference  material  exists.  Lengthy  mathematical  details  which 
are  of  minor  import  in  understanding  the  various  developments,  as  well 


as  certain  examples  and  error  analyses,  have  been  relegated  to  the 
appendices . 

In  Section  II,  the  exact  theory  describing  the  small-signal 
response  of  fluid  transmission  lines  of  annular  cross  section  is 
presented.  The  results  are  compared  with  those  obtained  by  others  for 
circular  and  rectangular  lines.  Thus,  Section  II  serves  also  as  an 
introduction  to  the  synthesis  of  results  applicable  to  circular  and 
noncircular  lines.  High  and  low  frequency  approximations  are  developed 
for  the  annular  line,  and  experimental  results  are  presented  which 
confirm  the  theory. 

In  Section  III,  a characteristic  frequency  and  corresponding 
characteristic  radius  are  determined,  which  allow  the  frequency  response 
of  noncircular  lines  to  be  presented  in  terms  of  equivalent  circular 
lines.  The  results  are  compared  with  those  previously  obtained  for 
rectangular  and  annular  lines  by  using  the  classical  hydraulic  radius 
together  with  circular  line  theory.  The  results  predicted  using  the 
characteristic  radius  and  circular  line  theory  are  shown  to  be  clearly 
superior  to  the  results  predicted  using  the  classical  hydraulic  radius 
and  circular  line  theory.  However,  the  hydraulic  radius  may  be  adjusted 
to  allow  its  use  in  the  circular  line  theory.  A concept  of  impedance 
matching  is  also  developed  in  Section  III. 

Many  practical  problems  are  more  concerned  with  the  time  domain 
response  of  transmission  lines.  In  Section  IV,  the  results  of  others, 
which  describe  the  impulse  and  unit  step  response  of  cii'culai'  lines 
onlv,  are  adapted  to  lines  of  noncircular  ci'oss  section  and  also 
compared  with  existing  experimental  data. 


Section  V is  devoted  to  the  simplification  and  synthesis  of 
results  describing  the  combination  of  steady  developed  turbulent  mean 
flow,  and  laminar  oscillatory  flow  produced  bv  a small  sinusoidal 
pressure  input  signal.  The  solutions  are  extended  to  include  lines 
of  noncircular  cross  section. 

In  Section  VI , the  frequency  domain  results  of  Section  V are 
transposed  into  the  time  detain,  and  the  various  combinations  of 
downstream  pressure  and  flow  response  to  upstream  step  inputs  of 
pressure  or  flow  are  presented. 

In  Section  VII,  the  last  section  of  this  work,  a summary  of  the 
conclusions  and  recommendations  is  presented. 
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1 1 . The  Small-Signal  Response  o:  Annular 
Fluid  Transmission  Lines 

Introduction 

The  small-signal  frequency  response  of  fluid  transmission  line 
of  circular  cross  section  has  been  obtained  in  various  forms  [1-3]*. 
Nichols  [3]  expressed  the  results  in  a manner  similar  to  that  used  for 
electrical  transmission  lines.  Schaedel  [k]  followed  closely  the 
approach  used  by  Nichols  and  obtained  the  response  for  lines  of 
rectangular  cross  section.  In  both  of  these  studies,  a distributed 
parameter  line  model  was  developed  and  results  were  obtained  for  the 
impedance,  admittance,  propagation  operator  and  characteristic  impedance 
of  the  line. 

There  are  a number  of  applications , such  as  aircraft  pitot-static 
systems,  hydraulic  mining,  and  instrumentation,  where  annular  lines  are 
used  and  the  small  signal  frequency  response  characteristics  are  of 
interest.  The  purpose  of  this  section  is  to  present  the  small -signal 
response  of  rigid  pneumatic  transmission  lines  of  annular  cross  section. 
Figure  1.  The  analytical  approach  taken  is  similar  to  that  used  by 
Nichols  [3]  and  Schaedel  [•*].  The  detailed  analytical  development  is 
presented  in  Appendix  C.  Finally,  the  annular  line  results  will  be 
compared  with  those  obtained  for  circular  and  rectangular  lines. 


lNumbers  in  brackets  designate  Bibliography  at  end  of  dissertation. 
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Frequency  Response 

Small  Signal  Equations  and  Solution.  The  equations  expressing 
conservation  of  mass,  momentum,  and  energy  for  laminar  small  signal 
oscillatory  flow  in  a pneumatic  line  of  annular  cross  section  are 
those  used  by  Nichols  [3]  to  obtain  a frequency  response  solution 
for  lines  of  circular  cross  section:  The  development  of  these 

equations  is  given  in  Appendix  C. 

3p  3(pu) 

3t  “ 3z 


v a r iii  - 1 3P 

r 3r  [_r  3rJ  + 3t  p 3z 

_^T  3_  f £T  1 _ 3T_  _ 6p  T 3P 
r 3r  3rJ  3t  p cD  3t 


The  dependent  variables  are  assumed  to  be  of  the  following  form 


T(r,z,t)  = T + T'  = T + AT(r,z)  e3 


u(r,z,t)  = u + u'  = u + u(r,z)  e' 


P ( z , t ) = P + P'  = P + AP( z) 


where,  for  example,  T is  the  average  temperature  and  T'  is  the  instanta- 
neous temperature  fluctuation.  Substitution  of  these  variables  into 
equations  (2)  and  (3)  gives 


3"u  1^  _3u  jeou  _ _1_  3AF 
gr2  r 3r  v “ vp  3z 


(8) 


3-at 

dr* 


i ^AT  juiAT  _ juiAP 

r 3r  v*|>  PCpV-j' 


Assuming  isothermal  walls  and  zero  velocity  at  the  walls,  the  boundary 
conditions  for  the  annular  line  are 

AT  = 0,  u = 0 at  r = r^  (9) 

AT  = 0,  u = 0 at  r = rQ  (10) 

The  momentum  and  energy  equations  are  thus  identical  in  form,  and 
the  solutions  for  the  dimensionless  velocity  and  temperature  profiles 
are  Bessel  functions  of  zero  order. 


The  arbitrary  coefficients  (B,D,E,  and  F)  obtained  from  the  boundary 
conditions  are 


(13) 
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The  coefficients  E and  E of  equation  (12)  are  identical  to  B and  D, 
respectively,  except  v is  replaced  by  v-j>  in  the  arguments  of  the 
Bessel  functions. 

Series  Impedance  of  the  Annular  Line.  The  mass  rate  of  flow  at 
the  reference  pressure  P*  is  given  by 


(pu)  2nrdr 


and  by  analogy  to  electrical  transmission  line  theory 


- = ZQ*  (1G) 

where  Z is  the  frequency  dependent  series  impedance  per  unit  length, 

AP  is  analogous  to  voltage  and  Q*  is  analogous  to  current. 

Solving  for  (pu)  from  equation  (11),  performing  the  integration  of 
equation  (15),  and  using  equation  (16),  an  expression  for  Z is  obtained: 


s.  & 


where 


M = B 


roJ 


IT ro^  • riJl(  ri^ 


jD 


r0Ki 
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(16) 


Introducing  the  adiabatic  inductance  per  unit  length,  La  = p*/An, 
equation  (17)  becomes. 


s jwL/ 


1 + 


2itM 


'»  $ 
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(19) 


where  Z is  shown  to  be  frequency  dependent. 

Shunt  Admittance  of  the  Annular  Line.  Trom  electrical  transmission 
line  theory. 


8Q* 

3z 


= YAP 


(20) 


where  Y is  the  shunt  admittance  per  unit  length.  Upon  differentiating 
equation  (15)  with  respect  to  z,  using  equation  (1),  and  substituting 
into  equation  (20),  Y becomes 


fXV 


Y = 


2t 

p*AP 


3p  , 

9t  rdr 


(21) 
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It  is  shown  in  Appendix  C that 


i£  B ££  II  + e p SL  . ill  (22) 

at  y at  p p I pep  at  at 

Introducing  the  complex  time  notation  for  3P/3t  and  3T/3t  and  using  the 
relation 

6pT 

CP-  cv  = -^  <23> 

equation  (22)  becomes 

3d  nr  r (y-1)  P Cp  AT  *| 

3t  = ^ j“AP  [X  + (Y-1} AP— j <**0 

where  6p  = 1/T  for  a perfect  gas.  Substituting  for  p Cp  AT/AP  from 
equation  (12),  equation  (24)  becomes 


Substituting  equation  (25)  into  equation  (21)  and  performing  the 
indicated  integration  yields 


where 


£P 

and  Ca  - ••  — ^ ■ - the  adiabatic  capacitance  per  unit  length.  For  an 
ideal  gas,  Ca  = An/yP*.  Tabulated  Bessel  functions  of  interest  are 
given  in  [5,6]. 

Equivalent  Electrical  Circuit.  The  series  impedance  and  shunt 
admittance  may  be  expressed  as 


= R + jioL 

(28) 

= G + 3 

(29) 

where  R,L,G,  and  C are  the  circuit  parameters  of  resistance,  inductance, 
conductance,  and  capacitance,  per  unit  length,  respectively.  The 
equivalent  circuit  for  a transmission  line  of  dz  length  is  shown  in 
Figure  2.  The  circuit  parameters  are  all  real  numbers  and  are  dependent 
on  the  signal  frequency,  the  cross-sectional  area,  and  the  radii  of  the 
annular  line  through  the  expressions  for  Z and  Y,  equations  (19)  and  (26). 

Characteristic  Impedance  and  Propagation  Operator.  Equations  (16) 
and  (20)  are  used  to  obtain  the  governing  differential  equations  for 
pressure  and  flow.  Differentiating  equation  (16)  with  respect  to  z and 
using  equation  (20)  gives 

ZY(dP)  (30) 
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The  general  solution  for  the  pressure  is 


AP  = Cje"1"2  + C2e  +rz 


(31) 


where  r,  the  propagation  operator,  is  defined  as 


r = i^Y  = a + j6 


(32) 


where  a is  the  attenuation  per  unit  length  and  6 is  the  phase  angle 
increase  per  unit  length.  The  solution  for  the  volumetric  flow  is 


Q * k (Cie'ri  ‘ C2'+rs' 


(33) 


where  Zo»  the  characteristic  impedance,  is  defined  as 

Z0  = - 1/Yq  (34) 

and  Yo  is  the  characteristic  admittance. 

The  characteristic  impedance  and  propagation  operator  are 
generally  complex  quantities  and  are  convenient  for  describing  the 
pressure  and  flow  in  the  transmission  line.  An  excellent  review  of 
the  significance  and  use  of  Zq  and  T in  various  representations  of 
fluid  lines  for  system  analysis  is  given  by  Goodson  and  Leonard  [7], 

The  phase  velocity  is 


c 


u 

B 


(35) 


where  X is  the  wavelength  and  6 increases  2v  radians  per  wavelength. 
The  phase  velocity  c is  in  general  not  equal  to  ca,  the  adiabatic 


15 


4 


The  propagation  velocity  ratio  or  wavelength  ratio  is 


The  neper  attenuation  in  Xa/2iv  distance  is 

a ca/u>  = a Xa/2n  (43) 

To  convert  the  attenuation  from  nepers  to  the  more  common  decibel, 
the  following  equations  are  used: 

(aXa)  db  = 40t  (logjpe)  a ca/u> 

* 54.  575  a ca/u>  db  per  Xfl  (44) 


(qX)  db  = 54.575  a/6  db  per  X (45) 

Characteristic  Frequencies.  Following  Nichols  [3]  it  is  convenient 
to  define  two  characteristic  frequencies  which  will  aid  in  the 
presentation  of  the  results  and  facilitate  comparison  of  the  annular 
line  results  with  those  obtained  for  circular  and  rectangular  lines: 


uv 


8nv 

^n 


s 


8v 

r0:(l-r*2 ) 


(46) 


8nv^  wv 


(47) 


where  r*  = r^/r0  is  defined  as  the  annular  radius  ratio.  The  viscous 
characteristic  frequency  wv  is  used  to  form  the  dimensionless  frequency 
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ratio  o)/(jjv  which  is  used  in  presenting  the  numerical  results.  Because 
the  characteristic  frequency  uv  is  defined  to  he  proportional  to  the 
ratio  v/A,  the  frequency  dependence  of  all  sices  of  annular  lines  of 
the  same  radius  ratio  are  represented  by  a single  curve  when  the 
results  are  plotted  against  <d/o)V.  This  is  analogous  to  the  frequency 
dependence  of  all  sizes  of  circular  lines  and  of  rectangular  lines  of 
the  same  aspect  ratio  when  plotted  against  u>/wv. 

Numerical  Results  and  Discussion 


The  numerical  results  were  calculated  on  a CDC  6600  computer  using 
the  standard  properties  of  air  at  80°F  given  in  Table  I.  Flots  are 
given  for  frequency  ratios  w/wv  from  C.l  to  1000  and  annular  radius 
ratios  of  0,  0.1,  0.6,  and  0.8. 


TABLE  I 

Properties  of  Air  at  80° F [3] 


F* 

14.696  psia 

P* 

1.1017^10"  lbf-sec*/in.  4 

V* 

0.02432  in. ‘/sec 

VT* 

0.03435  in. ‘/sec 

0“ 

0.7060  dimensionless 

Y 

1.4017  dimensionless 

ca 

13,674  in. /sec 

£c 

386.086  in./sec‘ 

The  influence  of  the  frequency  and  radius  ratio  on  annular  line 
performance  is  easily  seen  in  Figures  3 through  6.  The  modulus  of  the 
admittance  ratio  given  by  equation  (39)  is  plotted  in  Figure  3 and  the 
phase  angle  of  the  characteristic  admittance  is  plotted  in  Figure  4. 

It  can  be  seen  that  at  high  frequencies  the  magnitude  of  the  character- 
istic admittance  approaches  that  of  the  lossless  line,  while  the  phase 
angle  approaches  zero.  At  low  frequencies  the  phase  angle  approaches 
45  degrees.  The  propagation  velocity  ratio  c/ca,  or  wavelength  ratio, 
\/Aa,  equation  (42),  is  plotted  in  Figure  5.  For  a given  flow  area  and 
frequency,  the  propagation  velocity  decreases  with  increasing  radius 
ratio.  The  db  attenuation  per  line  wavelength,  equation  (45),  is 
plotted  in  Figure  6.  For  given  fluid  properties,  flow  area,  and 
frequency,  the  attenuation  per  line  wavelength  increases  with  increasing 
aspect  ratio. 

The  circuit  parameters  R,G,L,  and  C are  normalized  with  respect  to 
suitable  reference  values  and  the  results  plotted  in  Figures  7 through 
10. 

The  ratio  R/Rvn  is  plotted  in  Figure  7,  where  Kvn  the  viscous 
laminar  steady  flow  (DC)  resistance  of  the  annular  line,  and  is 
defined  in  terms  of  r*  [8]. 


(48) 


Equations  for  determining  the  DC  properties  of  circular,  annular,  and 
rectangular  lines  are  derived  in  Appendix  A.  As  the  frequency  u>  approache 
zero,  R approaches  Rvni  however,  RVn  is  a pood  approximation  for  R when 
< 1 for  all  r*.  for  the  larger  values  of  rft  shown  in  Figure  7, 


the  range  of  applicability  of  the  constant  resistance  model  is  extende 
to  higher  values  of  w/u>v. 


Figure  8 shows  the  frequency  dependence  of  L/La  where  La  is  the 
adiabatic  inductance.  The  tendency  of  the  frequency  dependent  induct- 
ance to  approach  its  adiabatic  value  at  high  frequencies  is  clearly  shown. 
For  low  frequency  ratios,  say  w/wv  < 1,  the  inductance  is  primarily 
dependent  on  r*  and  essentially  independent  of  u>/u>v.  The  value  of  L/La 
at  low  frequency  ratios  varies  between  approximately  1.2  and  1.33, 
depending  on  r*. 

The  ratios  G/w^  Ca  and  C/Ca  are  plotted  in  Figures  9 and  1C.  In 
Figure  9 the  conductance  is  shown  to  be  essentially  zero  when  u>/wv  < 1 
for  all  r*.  For  values  of  > 10  the  conductance  increases 

significantly.  This  increase  is  more  pronounced  at  large  values  of  r*. 
Figure  10  shows  the  general  thermodynamic  behavior  of  the  model: 


ca  ^n  n 
yP 


where  y is  the  adiabatic  exponent  (approximately  1.4  for  air  at  standard 
conditions)  and  n is  the  polytropic  exponent.  At  low  frequencies 
y/n  = 1.4;  thus,  n - 1 , which  is  the  isothermal  exponent.  At  high 
frequencies  y/n  = 1.0,  thus,  n = y,  the  adiabatic  exponent.  The  process 
is  therefore  isothermal  at  low  frequencies  and  adiabatic  at  high 
frequencies.  As  r*  increases,  the  transition  from  isothermal  to 
adiabatic  behavior  occurs  at  higher  values  of  w/wv. 
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Figure  6.  Attenuation  Per  Line  Wavelength  for  Annular  Lines 


Figure  10.  Ratio  of  Frequency  Dependent  to  Adiabatic  Capacitance 


In  summary,  at  low  frequencies  the  equivalent  circuit,  see 
Figure  2,  consists  essentially  of  constant  value  elements  of  resistance, 
inductance,  and  capacitance  and  no  conductance.  At  high  frequencies, 
the  inductance  and  capacitance  approach  their  adiabatic  values;  however, 
the  resistance  and  conductance  are  very  large  and  are  strongly  dependent 
on  frequency  and  radius  ratio. 

Comparison  with  Circular  and  Rectangular  Lines 

The  results  obtained  for  the  annular  line  are  now  compared  with 
those  obtained  for  circular  and  rectangular  lines.  The  comparison  is 
made  on  the  basis  that  the  respective  cross-sectional  areas  are  equal. 

Generally,  the  annular  results  approach  those  for  circular  lines 
as  the  inner  radius  of  the  annulus  approaches  zero  and  the  losses 
attributable  to  the  inner  wall  become  small.  It  can  be  shown  that  the 
circular  line  velocity  and  temperature  profiles  of  Nichols  [3]  may  be 
obtained  directly  from  equations  (11)  and  (12)  for  the  special  case  of 
the  annular  line  with  r*  = 0.  The  Bessel  function 


is  infinite  at  the  origin  and  ceases  to  be  a solution  in  equations  (11) 
and  (12).  Thus,  when  r^  = 0,  the  constants  D and  F are  identically 
zero.  Applying  the  boundary  conditions  of  equation  (10)  at  the  outer 
wall  r = rc  and  obtaining  new  relationships  for  B and  E leads  to 
Nichols'  expressions  for  the  velocity  and  temperature  profiles.  In 
all  of  the  figures  the  curves  for  r*  = 0 are  identical  with  those  of 
Nichols  for  circular  lines. 
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The  annular  results  converge  to  those  of  Schaedel  [4]  for  rectangular 
lines  as  r*  ■*  1.  The  comparison  is  made  by  defining  an  equivalent  aspect 


' 


ratio  for  the  annular  line.  The  aspect  ratio  of  the  rectangular  line  is 
defined  here  as  ARr  = b/h  2 1.  For  equal  rectangular  and  annular  flow 
areas, 

bh  = ir(r02-ri2)  = ir(r0+ri ) (r0-ri ) (50) 

where  choosing  b = iT(r0+ri)  and  h = r0-r^,  the  annular  aspect  ratio  is 
defined  as 


_ 11  ( ro+ri ) _ „(  1+r* ) 
n (rc-ri ) “ ( 1-r* ) 


(51) 


Since  the  annular  line  side  walls  of  dimension  h do  not  physically 
exist  and  are  not  considered  in  the  theory,  greater  losses  are  expected 
in  rectangular  lines  when  compared  to  annular  lines  of  the  same  area 
and  aspect  ratio.  However,  as  aspect  ratio  becomes  large,  the  losses 
attributable  to  the  short  walls  of  the  rectangular  line  become  small 
compared  to  the  losses  on  the  long  walls.  Consequently,  the  results  for 
annular  and  rectangular  lines  converge  with  increasing  aspect  ratio. 
Figure  11  is  a comparison  of  the  attenuation  per  unit  wavelength  (aX) 
of  annular  and  rectangular  lines  where  here  wv  is  defined  to  correspond 
to  that  used  by  Schaedel  [4];  i.e.,  uiv  = 4v/A.  Of  interest  is  the 
rather  close  agreement  at  moderately  low  aspect  ratios.  Agreement 
within  10%  is  achieved  for  all  frequency  dependent  properties  when 
AR  2 7 (r*  > 0.38).  Thus,  either  model  will  give  approximate  predictions 
of  the  performance  of  the  other  model  when  AR  > 7. 
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The  differences  in  the  circular,  rectangular,  and  annular  line 
results  for  equal  areas  are  due  mainly  to  differences  in  the  steady- 
flow  (DC)  resistances.  The  differences  in  the  DC  resistances  are 
illustrated  by  the  DC  resistance  ratios  plotted  in  Figure  12.  The 
viscous  resistance  per  unit  length  of  circular  lines  is  the  classical 
Hagen-Poiseuille  flow  resistance 


^vc 


8itu 


(52) 


For  equal  flow  areas  An  = Ac, 


^vn  _ 1+r*2  1 

^"vc  l-r*2  3-11  r*‘ 


(53) 


where  is  given  by  equation  (48).  The  ratios  Rvr/Rvc  311(3  Rvn/Rvr 
are  also  given,  where  is  the  viscous  resistance  of  rectangular  lines 
given  by  Schaedel  [4],  The  ratios  Rvn/Rvc  311(3  Rvr/Rvc  are  nearly 
linear  functions  of  AR  for  AR  >4.  The  convergence  of  Kvn/Kvr  to  unity 
for  large  AR  is  easily  seen. 

In  summary,  the  small  signal  performance  of  annular  lines  has  been 
presented  and  includes  the  performance  of  circular  lines  as  a special 
case  (r*  =0).  The  small  signal  performance  of  annular  and  rectangular 
lines  converge  with  increasing  aspect  or  radius  ratio.  Useful  approxi- 
mations of  rectangular  line  performance  may  be  made  using  the  annular 
line  results  at  higher  aspect  ratios,  or  vice  versa. 

It  is  important  to  note  that  whereas  the  numerical  results  are 
presented  for  air  at  80°F,  any  other  fluid  could  have  been  used.  If 


9 

\ 


s 
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the  fluid  is  considered  incompressible,  the  results  are  easily  adapted, 
and  considerably  simplified  by  assuming  y = 1,  and  by  taking  the  speed 
of  sound  as  c = »'kg/p  where  Kg  is  the  bulk  modulus  of  the  fluid.  It  is 
noted  that  the  shunt  admittance,  equation  (26),  becomes  a pure 
capacitance  of  constant  value  when  y = 1.  The  limitations  given  by 
Iberali  [1]  apply  to  the  annular  line  also: 

r0-r±  > — (=  1. 8*10“ 6 inches,  for  air) 
ca 

ca2 

v < (=  7.7*109  rad/sec,  for  air) 

v 

The  annular  line  model,  as  well  as  the  circular  and  rectangular  models 
of  [3]  and  [4],  respectively,  are  restricted  to  laminar  flow. 
Discrepancies  are  to  be  expected  for  turbulent  mean  flow  conditions. 

The  case  of  turbulent  mean  flow  is  considered  in  Section  V. 

High  and  Low  Frequency  Approximations 

Useful  high  and  low  frequency  approximations  of  the  system  variables 
are  most  readily  obtained  from  analysis  of  the  equivalent  circuit 
components  in  conjunction  with  the  applicable  Bessel  function 
approximations . 

For  the  low  frequency  approximations  (w/wv  <<  1),  the  equivalent 
circuit  consists  of  the  constant  value  elements  RVn>  Wn*  and  C^,  with 
Gi  = 0.  Here  is  the  isothermal  capacitance,  = yCa.  The  low 

frequency  series  impedance  and  shunt  admittance  are  respectively. 


ZLn  = Rvn  + 

YLn  = Gi  + j“>Ci  = ju>yCa 


( 54  ) 
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Using  the  inequality  wliVn/RVn  « 1 and  equations  (34)  and  (64),  the 
low  frequency  characteristic  admittance  and  propagation  operator  become 


Y0t  = ^'Ln/^Ln 


HsZ 

1 »vn 


<i+j) 


(55) 


rL  = ‘-Ln^Ln  = Q + <1+J>  (&6) 


It  is  seen  that  the  phase  lag  6 is  approximately  numerically  equal  to 
the  attenuation  a,  and  that  the  phase  angle  of  the  admittance  is 
approximately  r/4  radians.  For  u/wv  a 0.1,  equations  (55)  and  (56) 
give  F and  Yq  vithin  one  percent  of  the  exact  value,  for  all  r*. 
Expressions  for  Rvn  and  Lvr,  are  derived  in  Appendix  A. 

The  high  frequency  approximations  are  derived  in  Appendix  B using 

I 

the  standard  asymptotic  expansions  of  the  Bessel  functions  [5,6], 
along  with  other  approximations.  For  u/wv  s''  1,  the  inequalities, 

R/uiL  <<  1 and  G/wC  <<  1,  apply.  The  high  frequency  series  impedance 
and  shunt  admittance  become 

ZKn  = juLa„  j 1 . ^ 1 

'H„  - i“<4»  jl  ♦ O'r)  TP  (58) 

The  high  frequency  propagation  operator  and  characteristic  admittance 
become 
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rH 


aH  + 


1+r*  ' 
1-r* 


*U)WV  + j 


U) 


(59) 


It  is  of  interest  that  the  expression  for  in  equation  (59)  for 
annular  lines  is  identical  to  that  obtained  by  karam  and  Franke  [9]  for 
circular  lines  for  the  special  case  r*  = 0.  The  expression  for  is 
identical  to  that  given  in  [9].  It  should  therefore  be  possible  to 
obtain  a solution  in  the  time  domain  for  annular  lines  using  Karam's 
appi'oach  [19, PC].  The  frequency  to  time  domain  transformations  for 
circular,  rectangular  and  annular  lines  are  covered  in  Section  111. 

Experimental  Confirmation  of  Annular  Theory 

The  experimental  verification  of  the  annular  lino  theory  was  the 
primary  purpose  of  a M.S.  thesis  completed  by  Captain  G.  R.  Tarney 
at  the  Air  Force  Institute  of  Technology  [10].  Farnoy  measured  the 
small  signal  amplitude  and  phase  shift  frequency  response  of  a blocked 
annular  pneumatic  transmission  line  at  five  different  radius  ratios. 

The  outer  line  was  approximately  three  feet  long  and  had  an  internal 
diameter  of  0.187  inches.  The  five  radius  ratios  were  achieved  by 
inserting  and  centering  five  different  sired  rods  within  the  basic 
line,  rarney  observed  that  the  experimental  gains  wore  within  * 1 db 
and  the  phase  shift  angles  within  115  degrees  ol  those  predicted  bv 
the  annular  line  theory.  Further,  the  predicted  analytical  results  toll 


within  the  measurement  capability  of  the  experimental  instrument  at  ion  used . 


A description  of  the  experimental  apparatus  and  test  procedures  used 
is  given  in  [10]. 

Additional  equations  are  required  to  account  for  the  sending  and 
receiving  conditions.  These  are  standard  equations  of  electrical 
transmission  line  theory,  available  in  many  references  such  as  [11] 
and  are  repeated  herein  for  the  convenience  of  the  reader. 

For  a given  line,  the  transfer  pressure  gain  is  [12]: 


Zr  Zr*  . 

— cosh  T£  - ttt  sinh  Tl 
"s  "0 


(61) 


where  r and  s denote  sending  and  receiving  locations  and  i is  the 


distance  between  these  locations.  The  sending  impedance  is 

Zr  cosh  Ft  + Z0  sinh  n 
"s  = °°  Z0  cosh  + Zr  sinh  Tl 


(62) 


For  a blocked  line,  the  receiving  impedance  Zr  is  infinite,  thus 

Zs  = Z0  coth  (63) 

The  phase  lag  angle  between  sender  and  receiver  is 


6£ 


arctan 


Im_a£  (&)_ 
Real  (g) 


(6*4) 


The  pressure  gain,  equation  (61)  converted  to  decibels,  and  the  phase 
lag  angle,  equation  (6*+1  in  degrees,  are  plotted  as  functions  of 
excitation  frequency  in  Figures  13,  l*t,  and  15,  for  annular  radius 
ratios  of  0.07,  0.23,  and  0.67,  respectively.  The  experimental  data 
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rigure  1*4 • Correlation  of  Experimental  Results  with 
Annular  Line  Theory 


rigure  15.  Correlation  of  Experimental  Results  with 
Annular  Line  Theory 


of  Farney  is  also  plotted.  The  sharp  increase  in  phase  angle  shown 
in  the  figures  (approximately  720  degrees)  is  a plotting  artifice  to 
conserve  space  and  has  no  other  significance.  The  ambient  test 
conditions  at  each  radius  ratio  are  given  in  Table  II. 

TABLE  II 

Ambient  Test  Conditions'*’ 


It 

P( ambient) 

P(line) 

r 

(psia) 

(psig) 

0.07 

14.50 

5.01 

0.23 

14.44 

5.01 

0.67 

14.20 

5.01 

+The  tabulated  values  were  used  to  determine 
the  necessary  mean  properties  of  air. 

As  can  be  seen  from  the  figures,  the  agreement  between  theory  and 
experiment  is  excellent. 
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III.  Correlation  of  the  Linear  Response  of  Circular  and 
Noncircular  Fluid  Transmission  Lines 

In  steady  flow,  it  has  been  particularly  useful  to  correlate  the 
resistances  of  lines  of  noncircular  cross  sections  with  the  resistance 
of  circular  lines  through  the  classical  hydraulic  diameter.  Similarly, 
in  oscillatory  flow,  it  would  appear  useful  to  be  able  to  correlate  the 
frequency  response  of  noncircular  lines  to  the  response  of  circular 
lines.  Besides  the  comparisons  made  between  rectangular  and  circular 
lines  [4]  and  between  rectangular,  annular  and  circular  lines  herein, 
other  attempts  have  been  made  to  predict  the  response  of  rectangular 
lines  using  the  concept  of  the  hydraulic  diameter  [13,14], 

In  this  section,  a dimensionless  frequency  is  defined  which  leads 
to  a generalized  presentation  of  the  frequency  dependent  response  of 
lines  of  circular  and  noncircular  cross  section.  The  dimensionless 
frequency  in  turn  leads  to  a characteristic  radius  for  the  noncircular 
line  which  is  roughly  analogous  to  the  classical  hydraulic  radius  used 
in  steady  developed  flow  for  noncircular  lines.  It  is  shown  that  the 
frequency  dependence  of  noncircular  lines  may  be  very  closely  approxi- 
mated when  the  characteristic  radius  is  used  together  with  the  circular 
line  theory  of  Nichols  [3].  The  results  are  compared  with  those 
obtained  by  Healey  and  Carlson  [13]  for  rectangular  lines.  In  [13],  the 
classical  hydraulic  radius  together  with  Nichols’  theory  was  used  to 
predict  the  response  of  rectangular  lines.  The  response  predicted 
using  the  characteristic  radius  defined  herein  is  shown  to  be  clearly 
superior  to  the  response  predicted  using  the  hydraulic  radius.  Results 
are  presented  for  both  rectangular  and  annular  lines  over  a wide  range 


42 


— 


of  aspect  and  radius  ratio,  respectively.  In  addition,  an  approach  to 
impedance  matching  based  on  equivalent  circular  lines  is  presented. 

Frequency  Response  Solutions 

Several  dimensionless  frequencies  (or  frequency  ratios)  have  been 
defined  in  order  to  present  generalized  frequency  response  solutions 
for  circular  lines.  Nichols  [33  showed  that  the  frequency  response  of 
a circular  line  is  a function  of  the  dimensionless  frequency  ratio 
w/uv,  where  u;  is  the  angular  excitation  frequency  and  wv  is  a viscous 
characteristic  frequency  defined  as 


Ry 

L* 


(65) 


Rv  is  the  laminar  DC  viscous  resistance  per  unit  line  length  and  La  is 
the  adiabatic  inertance  per  unit  line  length.  For  any  line,  Rv  and  La 
are  defined  as 


(66) 


and 


La  = 7T 


(67) 


For  a circular  line,  Rv  and  u>v  become 


R - 

Rvc  " . 2 


(68) 


and 


U)v 


8nv 


C hr 


(69) 


t 
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Other  formulations  of  the  dimensionless  frequency  have  also  been 
defined  for  the  circular  line;  for  example,  Brown  et  al.  [15]  used 
fl  = tor2/v  = 8u)/o)nic,  while  Goldschmied  [16]  used  SN  = 4usr2/\>  = 32  w/u>Vc> 
In  all  of  these  cases,  however,  the  dimensionless  frequencies  are 
proportional  to  toA/v. 

Schaedel  [4]  showed  that  the  response  of  rectangular  lines  of  the 
same  aspect  ratio  is  a function  of  tuA/4 v,  and  it  has  been  shown  herein 
that  the  response  of  concentric  annular  lines  of  the  same  radius  ratio 
is  a function  of  wA/8mv.  Typical  results  for  the  attenuation  per  line 
wavelength  of  rectangular  and  annular  lines  are  plotted  in  Figures  16 
and  17,  respectively.  Circular  line  results  are  included  in  Figure  17 
as  a special  case  of  the  annular  line  with  the  inside  radius  r^  = 0. 

The  dependence  on  aspect  ratio  and  radius  ratio  as  well  as  signal 
frequency  is  clearly  shown. 

Generalized  Response  of  Non circular  Lines 

The  a A curves  shown  in  Figures  16  and  17  for  different  aspect 
ratios  and  radius  ratios  all  have  the  same  general  shape.  Thus,  it 
would  seem  possible  that  the  results  for  different  aspect  or  radius 
ratios  could  be  generalized  and  related  to  a single  curve.  This  can  be 
accomplished  approximately  by  using  the  dimensionless  frequency  defined 
by  equation  (65)  with  Rv  taken  as  the  resistance  of  the  particular 
cross  section  being  considered. 

A better  approach  is  to  include  the  inertance  effects  in  the 
dimensionless  frequency  as  well  as  the  resistance  effects  and  to  present 
the  results  as  a function  of  the  dimensionless  frequency  u>/wc,  where  u>c 
is  the  radian  frequency  defined  by  Nichols  [3]  as 


(70) 


Here  Rv  is  defined  by  equation  (66),  and  Lv,  the  low-frequencv  or 
viscous  inertance,  is  defined  by 


Lv 


dA 


La  = 


(71) 


where  depends  on  the  laminar  fully-developed  velocity  distribution. 

The  ratio  oj/wc  includes  both  viscous  and  inertial  effects  at  low 
frequency.  Consequently,  this  ratio  more  accurately  represents  the 
fluid  line  at  low  frequency  ratios,  and  gives  analytical  results  which 
are  in  closer  agreement  with  both  the  theoretical  and  experimental 
results  obtained  for  circular  and  noncircular  lines. 

The  a\  results  shown  in  Figures  16  and  17  are  replotted  in 
Figure  18  as  a function  of  u/wc.  The  families  of  curves  for  different 
aspect  and  radius  ratios  are  now  reduced  approximately  to  a single 
curve.  By  incorporating  the  geometric  dependence  of  aspect  and  radius 
ratio  in  wc,  lines  of  different  sines  and  cross  sections  are  shown  to 
have  approximately  the  same  generalined  frequency  behavior.  The 
circular  line  attenuation  is  also  included  in  Figure  18,  where  wc  for 
the  circular  line  is  the  same  as  that  given  by  Nichols  [3D;  i.e., 

8nu 


u>Cc  - 


ac 


61TV  3 

Ac  = 4 Wvc 


(72) 


Figure  18  shows  aX,  the  attenuation  per  line  wavelength.  In  order  to 
obtain  a,  the  attenuation  per  unit  length,  it  is  necessary  to  know  X 
the  wavelength  in  the  line  for  each  ratio  w/wc.  Curves  of  X/Xa  vs 
u>A/4v  for  rectangular  lines  are  given  in  [4],  and  are  not  reproduced 
herein.  However,  for  annular  lines  curves  of  X/Xa  vs  uxA/8*v  are  given 
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Attenuation  Per  Line  Wavelength  for  Circular,  Rectangular,  and 


I 


? 


p* 


in  Figure  5.  The  dependence  on  radius  ratio  is  clearly  seen.  Similar 
qualitative  dependence  on  aspect  ratio  may  be  seen  in  [4]  for  rectangular 
lines.  The  families  of  curves  depicting  X/Xa  reduce  approximately  to 
one  curve  representing  circular,  rectangular  and  annular  lines  wher 
plotted  against  io/u)c  as  shown  in  Figure  19.  The  ratio  X/Xa,  equation 
(42),  is  also  the  ratio  of  the  propagation  or  phase  velocity  in  the  line 
to  the  free  sonic  (adiabatic)  velocity;  i.e.,  X/Xa  = c/ca. 

Numerical  Results  and  Discussion 

The  results  in  Figures  16  through  19  were  computed  assuming 
standard  properties  of  air  at  80  F (Table  I).  The  numerical  differences 
in  aX  in  Figure  18  for  the  range  of  aspect  or  radius  ratios  shown  in 
Figures  16  and  17  were  small  and  are  not  plotted  in  Figure  18.  The 
maximum  difference  in  aX  between  circular  and  rectangular  lines  was 
1.05  db  at  w/u>c  = 6.0  for  a = 10.  The  maximum  difference  between 
annular  and  circular  lines  was  1.40  db  at  u>/wc  = 6.0  for  r“  = 0.8.  For 
u)/toc  < 0.1  the  difference  was  less  than  0.09  db,  and  for  w/wc  > 10  the 
difference  was  less  than  1.25  db. 

The  maximum  difference  in  X/Xa  between  circular  and  rectangular 
lines  was  C.035  at  u/( ioc  = 2.0  and  a = 10.  The  maximum  difference 
between  circular  and  annular  lines  was  0.046  at  w/u)c  = 3.0  and  rA  = 0.8. 
For  u)/wc  < 0.1  the  maximum  difference  for  both  rectangular  and  annular 
lines  was  less  than  0.017  and  for  id/u)c  > 10  the  maximum  difference  was 
less  than  0.025.  The  maximum  departure  of  X/Xa  from  that  predicted  for 
the  circular  line  for  the  range  of  aspect  and  radius  ratios  considered 
occurred  for  the  annular  line  of  radius  ratio  r*  = 0.8,  as  shown  in 
Figure  19.  All  other  configurations  are  nearer  to  the  circular  line 
results . 
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Figure  I'i.  Wavelength  catio  for  Circular,  Pectangular,  and  Annular  Lines 


It  is  noted  that  whereas  aX  vs  w/wc  gives  the  tightest  grouping  of 
the  aX  curves,  X/Xa  vs  w/wv  gives  a slightly  better  grouping  of  the 
A/Xa  curves,  where  wv  is  defined  by  equation  (65). 

Determination  of  the  Characteristic  Frequency 

The  determination  of  values  of  tuc  for  noncircular  lines  reduces 
essentially  to  obtaining  values  of  Rv  and  Ly.  Analytical  expressions 
of  Rv  and  Lv  for  circular,  rectangular  and  annular  lines  are  developed 
and  discussed  in  Appendix  A.  A convenient  method  for  obtaining  Rv  is 
to  relate  the  Rv  of  the  line  with  the  noncircular  cross  section  to  that 
of  a circular  cross  section;  i.e.. 


where  Kr  for  rectangular  lines  and  annular  lines  depends  on  aspect 
ratio  and  radius  ratio,  respectively.  Values  of  KR  have  been  calculated 
[4,17]  and  a few  are  given  in  Table  III.  The  viscous  inertance  Ly  is 
obtained  from  values  of  Kl»  equation  (71),  for  each  geometric 
configuration.  Values  of  have  also  been  calculated  [4,17,18]  and  a 
few  are  given  in  Table  III. 

Starting  with  equation  (70),  a convenient  equation  for  computing 
wc  is  easily  derived  as 


Kr  8trv 
“c  = KL  A 


(74) 


For  circular  lines:  Krc  = 1 and  Klc  = 4/3;  thus  wCc  = 6ttv/Ac.  Values 
of  wc  for  noncircular  lines  can  easily  be  obtained  using  values  from 
Table  III  or  using  the  equations  or  data  from  the  references  [4,17,18]. 
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TABLE  III 


Fully  Developed  Laminar  Flow  Parameters  for 
Rectangular  and  Annular  Lines 


KRr 

Kr 

Lr 

frRer 

ri/rO 

K*n 

Kr 

^n 

^n^en 

1.1322 

1.3785 

56.91 

1.2733 

0 

1.0000 

1.3333 

64.00 

1.0000 

1.3920 

1.3475 

62.19 

1.4324 

0.1 

1.7068 

1.2298 

89.37 

1.2222 

1.8133 

1.3120 

68.36 

1.6977 

0.2 

2.1645 

1.2168 

92.35 

1.5000 

2.2671 

1.2876 

72.93 

1.9895 

0.3 

2.7232 

1.2101 

93.85 

1.8571 

2.7317 

1.2714 

76.28 

2.2918 

0.5 

4.4649 

1.2035 

95.25 

3.0000 

5.0959 

1.2365 

84.68 

3. 8515 

0.8 

13.4890 

1.2004 

95.92 

9.0000 

oo 

1.2000 

96.00 

CO 

1.0 

oo 

1.2000 

96.00 

OO 

The  steady  flow  laminar  resistance  Rv  does  not  depend  on  frequency; 
consequently,  Rv  may  be  expressed  in  terms  of  the  classical  steady- 
flow  parameters: 


Rv  — fRe 


2Dh2A 


where  is  the  conventional  hydraulic  diameter,  f is  the  friction 
factor  and  Re  is  the  Reynolds  number.  For  a circular  line,  D^  = Dc 
and  fRe  = 64;  therefore, 


2Dc2AC 


52 


Dividing  equation  (75)  by  equation  (76) 


Rv_ 

Rv, 


fRe  Dc2  Ac 
64  A 


(77) 


where  fRe,  D^,  and  A are  determined  for  the  noncircular  cross  section 
being  considered.  Selected  values  of  fRe  are  given  in  Tab_e  III  for 
rectangular  and  annular  cross  sections.  The  relationship  between  Kr 
and  fRe  is  obtained  from  equations  (73)  and  (77) 


Kr  = 


fRe  ^c2  A_ 
64  ^2  Ac 


64  ^ 


(78) 


where  Kq  is  a geometric  shape  factor  defined  as 


Rg 


A_ 

Dh2  Ac 


4A 

*Dh2 


P 2 
rw 

4ttA 


(79) 


Kq  is  a constant  for  each  geometric  cross  section.  It  is  important  to 
note  that  when  = Dc  the  areas  of  the  noncircular  cross  section  A and 
equivalent  diameter  circle  Ac  are  not  equal  and  Kg  = A/Ac.  Likewise, 
when  A = Ac  the  hydraulic  diameters  are  not  equal  and  Kq  = Dc2/Dy)2. 
Expressions  for  Kq  are  given  in  Table  IV  for  circular,  rectangular  and 
annular  cross  sections.  Selected  values  of  Kq  are  included  in  Table  III 
for  rectangular  and  annular  cross  sections. 
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TABLE  IV 

Geometric  Cross  Section  Factor 


Cross  Section 

Kg 

circular 

i 

rectangular 

(a+l)2/ira 

concentric  annular 

l+r*/l-r* 

Use  of  the  Hydraulic  Diameter 

Results  based  on  hydraulic  diameter  can  be  used  to  predict  line 
performance  for  small  signal  oscillatory  flow.  At  a given  signal 
frequency.  Figure  18  shows  that  for  the  same  attenuation  on  lines  of 
different  cross  section,  the  <oc  must  be  the  same.  The  use  of  equivalent 
hydraulic  diameter  for  a noncircular  line,  however,  leads  to  an  wc  that 
is,  in  general,  not  equal  to  u>cc.  This  will  now  be  illustrated  for  a 
rectangular  line. 

From  equation  (74),  the  ratio  w^/wcc 
wcp  = 4 KRr  Ac 

“)cc  3 KLr  Ar  (80) 


By  substituting  for  Kj^r  from  equation  (78),  the  ratio  becomes 


^-r  _ 4 frRer  1 Dc2 
U)cc  3 64  ^Lp 


(81) 


For  Dhr  = Dc,  equation  (81)  reduces  to 
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(82) 


* 3K^~6lT 


Typical  values  of  <*)Cr/t*>Cc  and  *4»>Cr,/u)Ccl  as  a function  of  aspect 
ratio  are  shown  in  Table  V using  the  values  of  frRer  and  KLr  from 
Table  III. 


TABLE  V 

Frequency  Ratio  Parameters  for  Rectangular 
Lines  for  Dhr'  = Dr. 


b/h 

WCr/<*>CC 

•^Cp/UCc' 

1 

0.  86 

0.93 

2 

0.96 

0.98 

3 

1.09 

1.04 

4 

1.18 

1.09 

5 

1.25 

1.12 

10 

1.43 

1.19 

The  ratio  ^Cv/^cc  can  be  interpreted  approximately  as  the  ratio  of 
the  viscous  attenuation  of  rectangular  to  circular  lines  of  the  same 
hydraulic  diameter  at  a given  frequency.  The  values  of  /wCr/wCc  given 
in  Table  V confirm  the  attenuation  results  of  Healey  and  Carlson  [13] 
and  explain  the  differences  they  found  between  the  rectangular  line 
solutions  and  circular  line  solutions  of  the  same  hydraulic  diameter. 
When  using  circular  line  equations  to  predict  rectangular  line  results, 
the  differences  can  be  avoided  if  the  rectangular  line  friction  factor 
(frRer)  and  inertance  Lvr  are  included  in  the  calculations.  This  point 
will  be  made  clearer  in  the  next  few  paragraphs. 
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Determination  of  the  Characteristic  Radius 

It  is  important  to  note  that  uic  characterizes  the  frequency 
dependent  line  performance,  and  that  to  compute  wc,  only  the  kinematic 
viscosity,  line  geometry,  and  steady  flow  laminar  velocity  profile  need 
be  known.  Presuming  that  uic  = uicx  is  known  for  a line  of  arbitrary  but 
constant  axial  cross  section  an  equivalent  circular  line  containing  the 
same  fluid  is  sought.  The  characteristic  frequency  of  the  circular 
line  is 


u)Cc 


Rvc  _ 6v_ 
^c  ' r02 


(83) 


Setting  wCc  = wCx,  the  radius  of  an  equivalent  circular  line,  called 
the  characteristic  radius,  is  easily  found  to  be 

rc  = *^6v/ij)Cx  (84) 

Since  no  additional  restrictions  have  been  placed  on  the  arbitrary'  line 
(i.e.,  neither  the  flow  areas  nor  hydraulic  diameters  are  assumed  equal), 
equation  (84)  implies  that  for  each  arbitrary  line  there  exists  a unique 
equivalent  circular  line.  The  equivalent  circular  line  will  have 
approximately  the  same  normalized  frequency  dependent  behavior  as  the 
arbitrary  line  since  the  characteristic  frequencies  are  equal.  The 
characteristic  radius  is  roughly  analogous  to  the  hydraulic  radius  used 
in  steady  flow  problems. 
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Comparison  of  Characteristic  and  Hydraulic  Radii 

Healey  and  Carlson  [13]  used  the  classical  hydraulic  radius  (based 
on  the  hydraulic  diameter)  for  rectangular  lines. 


2A  _ Jbh  _ £h 
^ = Pw  ‘ b+h  2 


(85) 


together  with  Nichols'  equations  for  circular  lines  [33  to  estimate  the 
frequency  response  of  blocked  pneumatic  lines  of  rectangular  cross 
section.  These  results  were  then  compared  with  the  results  predicted 
by  the  rectangular  line  theory  of  Schaedel  [4]  and  with  their  own 
experimental  data.  The  experimental  data  and  predictions  using 
Schaedel' s theory  were  in  excellent  agreement  for  all  aspect  ratios 
tested  (a=l,2,3,4,5).  Using  rjj  and  Nichols'  equations  for  circular 
lines  yielded  good  agreement  with  the  data  and  with  Schaedel' s theory 
for  aspect  ratios  of  2,  3,  and  4;  however,  at  a=l  and  a=5  significant 
differences  were  observed  between  the  predicted  downstream  amplitude 
response  using  r^  and  the  measured  response. 

Using  the  characteristic  radius,  equation  (84),  together  with 
Nichols'  theory,  the  normalized  attenuation  aX  vs  oo  was  obtained  for 
aspect  ratios  a = 1,  3,  5,  and  10  and  the  results  plotted  in  Figure  20. 
For  comparison,  predictions  using  r^  and  Nichols'  theory  as  in  [13], 
and  using  Schaedel 's  theory  are  also  plotted.  The  differences 
between  Schaedel's  results  and  those  using  rc  are  so  small  that  they 
cannot  be  seen  on  Figure  20.  It  may  be  concluded  that  the  use  of  rc 
gives  much  better  agreement  with  Schaedel's  theory  and  consequently 
with  the  experimental  data  of  [13],  than  does  the  use  of  rj,.  It  can 
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igare  20.  Comparison  of  Attenuation  Per  Line  Wavelength  of  Rectangular  Lines 


also  be  seen  that  the  use  of  hydraulic  radius  predicts  toe  much  attenuation 
per  line  wavelength  for  a=l  and  not  enough  for  ab3.  The  difference 
increases  with  aspect  ratio  above  a=3. 

A similar  comparison  is  presented  for  annular  lines  of  r*  = 0, 

0.1,  0.5,  and  0.8  in  Figure  21.  At  r*=0,  rc  = r^  = r0  and  Nichols'  curve 
for  circular  lines  is  recovered.  For  r*>0,  the  attenuation  predicted 
using  and  Nichols'  theory  is  less  than  that  predicted  using  either 
rc  anc  Nichols'  theory’  or  the  exact  annular  theory  given  in  Section  II 
while  the  latter  two  approaches  are  in  excellent  agreement.  The 
difference  in  predicted  performance  based  on  rc  and  r^  continue  to 
diverge  with  increasing  radius  ratio. 

It  is  of  interest  to  compare  the  hydraulic  radius  and  the 
characteristic  radius  rc  of  a rectangular  line,  for  a practical  range 
of  aspect  ratio,  Figure  22.  The  two  approaches  yield  the  same 
equivalent  circular  line  at  a=2.3  (rc=r^,),  and  are  nearly  equal  over 
the  range  2 < a s 4.  However,  the  divergence  of  rc  and  r^  away  from 
a-2.3  is  clearly  seen.  The  circular  radius  required  to  give  the  same 
flow  area,  denoted  rA  is  defined  as 


and  is  used  to  normalize  both  r^  and  rc.  Similar  results  are  plotted 
for  annular  lines  in  Figure  23.  It  can  be  seen  that  the  ratio  rc/rh 
decreases  rapidly  from  unity,  where  r*=0,  and  rc  = r^  = rA.  The  area 
of  the  equivalent  circular  line  using  the  characteristic  radius  is 
approximately  60%  of  that  obtained  using  the  hydraulic  radius  for 
r^  > 0.2. 
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ius  Ratios 


In  principle,  if  the  frequency  dependence  of  lines  of  arbitrary 
axial  cross  section  may  be  correlated  as  above  using  uir,  then  an 
approximate  impedance  match  may  be  found  when  it  is  necessary  to  join 
lines  of  different  axial  cross  sections,  as  discussed  in  the  next  few 
paragraphs . 

Impedance  Matching 

It  is  of  practical  importance  to  terminate  a transmission  line  in 
its  characteristic  impedance  or  admittance  since  such  a termination 
implies  no  reflections.  Minimizing  reflected  signal  strength  is 
synonymous  with  increasing  the  efficiency  of  the  line  in  transmitting 
the  desired  signal.  Thus,  it  is  appropriate  to  discuss  impedance 
matching  when  lines  of  different  axial  cross  section  are  to  be  joined. 
Such  circumstances  occur  quite  often  in  fluidic  applications  where 
circular  lines  feed  fluidic  components  whose  basic  flow  regions  are 
rectangular  due  to  laminate  techniques  in  component  fabrication.  The 
joining  of  flow  regions  of  different  axial  cross  section  may  also  be 
found  in  aircraft  hydraulic  systems. 

If,  for  the  present,  we  ignore  the  junction  itself  and  concentrate 
on  matching  the  characteristic  impedance  of  transmission  lines  of 
different  cross  section,  it  would  seem  appropriate  to  require  the 
characteristic  frequencies  to  be  matched. 

The  modulus  and  phase  angle  of  the  characteristic  admittance  ratio 
for  circular,  rectangular  and  annular  lines  versus  u)/ucare  plotted  in 
Figures  24  and  25,  respectively.  Again  it  is  seen  that  lines  with  the 
same  wc  exhibit  approximately  the  same  frequency  dependence.  The 
admittance  ratio  of  rectangular  lines  is  within  4%  of  that  of  circular 
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Figure  21*.  Characteristic  Admittance  Patio  for  Circular,  Pectangular  and  Annular  Lines 
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lines  at  aspect  ratios  a < 10,  and  the  admittance  ratio  of  annular 
lines  is  within  6%  of  that  of  circular  lines  for  radius  ratios  r*  < 0.8. 
At  lower  aspect  or  radius  ratios,  these  percentages  decrease.  The 
phase  angle  of  the  admittance  vs  u>/ujc  for  both  rectangular  and  annular 
lines  is  within  1.5  degrees  of  that  of  circular  lines  for  the  same 
ranges  of  aspect  and  radius  ratios. 

It  is  important  to  note  that  the  modulus  of  the  characteristic 
admittance  has  been  normalized  with  respect  to  that  of  a lossless  line 
of  the  same  cross  section.  In  general  terms 


Yo 


‘sO 


“ f(oj/uc) 


or  |Y0|  a |Ys0|*f(«/»c) 


(87) 


and  equations  (87)  are  valid  for  all  lines.  It  follows  that  if  two 
lines  have  the  same  characteristic  frequency,  then  for  a given  excita- 
tion frequency  w the  lines  will  have  approximately  the  same  admittance 
ratio  but  not  necessarily  the  same  admittance.  From  the  definition  of 

Ysq  * 


Y 


so  = 


(36) 


and  Yso  is  obviously  area  dependent.  Therefore,  setting  the  character- 
istic frequencies  equal  does  not  insure  that  the  characteristic 
admittances  are  equal.  If  the  conditions  that  Yq  and  u)c  must  be  the 
same  are  imposed  on  two  lines  carrying  the  same  fluid,  then  the  flow 
areas  must  also  be  the  same.  These  conditions  are  in  general  incompat- 
ible except  for  the  special  case  when  the  two  lines  are  identical.  The 


conditions  amount  to  requiring  that  r^  = rc  when  lines  of  circular  and 
noncircular  cross  section  are  to  be  joined.  From  Figure  23,  it  is  seen 
that  for  annular  lines  r^  = rc  only  when  r*  = 0\  for  rectangular  lines. 
Figure  22,  this  condition  is  never  achieved.  It  is  concluded  that  the 
characteristic  admittance  (or  impedance)  of  lines  of  different  cross 
section  carrying  the  same  fluid  cannot  be  matched  at  each  and  every 
physical  frequency  ui.  The  closest  match  to  a circular  line  is  obtained 
when  r^=vc.  (rrt  + 0 or  a 4 1,  for  annular  or  rectangular  lines.) 

In  practice,  the  most  common  constraint  encountered  is  the 
specification  of  a steady  mean  flow  through  which  the  signals  are  to 
be  transmitted.  This  constraint  is  equivalent  to  requiring  that  the 
flow  continuity  is  maintained.  It  follows  that  for  a given  fluid  two 
lines  of  the  same  flow  area  have  the  same  lossless  admittance.  At  a 
specified  frequency  ratio  oj/<dc,  the  two  lines  will  also  have  approxi- 
mately the  same  frequency  dependent  characteristic  admittance.  However, 
the  characteristic  frequencies  are  not  the  same  and,  consequently, 
neither  are  the  physical  frequencies  for  which  the  admittances  are  equal. 

The  characteristic  frequency  of  noncircular  lines  will  be  higher 
than  that  of  circular  lines  of  the  same  cross-sectional  area.  For  any 
excitation  frequency  oi,  this  implies  that  the  characteristic  admittance 
(impedance)  of  noncircular  lines  will  be  lower  (higher)  than  for 
circular  lines  of  the  same  area,  and  the  designer  is  unfortunately 
faced  with  the  resultant  mismatch.  Nevertheless,  the  characteristic 
frequency  (or  radius)  of  the  noncircular  line  together  with  the  circular 
line  theory  may  be  used  in  a simple  procedure  to  predict  the  performance 
of  the  noncircular  line  at  any  desired  frequency. 
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The  ratios  u>cr/wCc  and  wCn/uiCc  are  plotted  in  Figures  26  and  27 
for  rectangular  and  annular  lines,  respectively,  assuming  equal  flow 
areas.  Similar  curves  could  be  constructed  for  any  noncircular  line. 
The  characteristic  frequency  of  circular  lines  is  easily  determined 
(«cc  = 6nv/Ac)  and  from  curves  similar  to  Figures  26  and  27  the  chara- 
cteristic frequency  (uc)  of  the  particular  noncircular  line  is 
determined.  For  the  physical  frequency  of  interest,  the  frequency 
ratio  u)/wc  is  computed  and  any  desired  performance  parameter  such  as 
aX,  Figure  19,  or  | Y q /Ys 0 1 , Figure  24,  is  easily  obtained  from  the 
circular  line  performance  curves.  It  is  emphasized  that  laminar  mean 
flow  is  assumed;  combined  turbulent  mean  flow  and  laminar  oscillatory 
flow  is  discussed  in  Section  V. 


IV.  The  Transient  Response  of  Noncircular  Fluid 
Transmission  Lines 

The  transient  response  of  fluid  transmission  lines  of  circular 
cross  section  was  obtained  by  Brown  £2},  who  derived  the  propagation 
operator  and  characteristic  impedance  of  an  infinitely  long  uniform 
rigid  line  in  the  Laplace  domain  for  both  compressible  and  incompress- 
ible fluids.  Brown's  solution  is  identical  to  that  of  Iberall  Cl]  and 
Nichols  [3]  frequency  domain  models  for  the  special  case  of  sinusoidal 
excitation  of  the  circular  line.  Subsequently,  Karam  [19,20]  developed 
a simplified  solution  for  the  transient  response  of  circular  lines 
after  transforming  approximate  high  and  low  frequency  domain  models 
into  the  time  domain  and  then  observing  a basic  similarity  between  the 
transformed  solutions.  Karam' s solution  is  more  easily  applied  than 
that  of  Brown  [2],  and  may  be  used  to  obtain  the  impulse  and  step 
responses  of  semi-infinite  fluid-filled  circular  lines  with  good 
accuracy  and  without  the  aid  of  a computer. 

In  Section  III,  a method  of  expressing  the  frequency  dependence  of 
noncircular  lines  entirely  in  terms  0f  equivalent  circular  lines  was 
presented,  and  results  were  obtained  for  rectangular  and  annular  lines 
over  a wide  range  of  aspect  and  radius  ratio.  The  primary  purpose  of 
this  section  is  to  present  a method  of  transforming  the  frequency 
response  of  noncircular  lines  into  the  time  domain.  The  analytical 
approach  essentially  combines  the  approaches  of  [2,3,20]  and  that  of 
the  previous  section,  and  results  in  a simple  but  complete  solution  of 
the  step  response  of  noncircular  lines. 
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frequency  Domain  Equations 


The  exact  analytical  expressions  for  the  parameters  Z,  Y,  r,  and 
Y0  which  characterize  the  performance  of  a fluid  transmission  line  are 
relatively  cumbersome  to  use  in  the  frequency  or  Laplace  domain. 
Furthermore,  it  has  been  observed  that  the  time  transforms  associated 
with  these  parameters  have  not  yielded  closed  form  analytical  expressions 
over  any  of  the  frequency  or  Laplace  domain,  even  for  the  simple  case 
of  circular  transmission  lines.  The  mathematical  difficulties  may  be 
traced  directly  to  the  underlying  frequency  or  Laplace  domain  solutions, 
which  contain  ratios  of  Bessel  functions.  These  ratios  are  not  thus 
far  expressible  in  a form  amenable  to  transformation  into  the  time 
domain.  Further,  according  to  Sneddon  [21],  it  is  unlikely  that  the 
Bessel  solutions  will  ever  be  totally  transformed.  Instead,  high  and 
low  frequency  approximations  are  made,  and  these  approximations  are  then 
transformed  into  the  time  domain.  In  the  mid- frequency  range,  very 
complicated  numerical  techniques  have  been  used  to  accomplish  the 
transformation  [22].  Alternately,  the  transient  response  in  the  mid- 
frequency range  may  be  approximated  in  the  time  domain  itself,  as  was 
done  by  Karam  [20]. 

The  general  approach  taken  herein  will  be  to  partition  the 
frequency  domain  into  three  regions , to  analyze  the  equivalent 
electrical  circuit  of  a general  transmission  line  in  each  of  these 
regions,  and  then  to  express  the  results  in  the  time  domain.  In  this 
approach,  the  text  by  Weber  [11]  is  especially  useful. 

The  three  frequency  regions  are  defined  with  reference  to  the 
characteristic  frequency,  o>c  = Rv/Lv*  as  follows:  low  frequency 
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(to  « wc),  mid- frequency  (to  = u)c),  and  high  frequency  (to  >>  toc).  The 
characteristic  frequency  ratio  to/toc  = toLv/Rv  may  be  thought  of  as  a 
frequency  dependent  (AC)  Reynolds  number.  At  low  frequency  ratios 
viscous  flow  effects  are  prevalent  while  at  high  frequency  ratios 
inertial  effects  prevail.  In  the  mid-frequency  range  viscous  and 
inertial  effects  are  of  the  same  order  of  magnitude.  Similarly,  the 
thermal  characteristic  frequency  lo-p  = toc/c2  roughly  divides  the  heat 
transfer  effects  into  regions  of  isothermal  behavior,  to  <<  and 
adiabatic  behavior,  to  » toq1.  Since  the  Prandtl  number  (o2)  is  near 
unity  for  air  (and  many  other  gases)  either  toc  or  uj  will  suffice  as 
the  overall  behavioral  index,  and  <oc  has  been  chosen  for  this  work. 

For  liquids  (y  = 1),  there  is  essentially  no  dependence  on  the  Prandtl 
number,  since  terms  in  the  energy  equation  containing  the  Prandtl 
number  are  premultiplied  by  the  quantity  (1-y).  [See,  for  example, 
equation  (26).] 

At  low  frequencies,  the  equivalent  circuit.  Figure  2,  is  comprised 
of  essentially  constant  value  elements  of  R,  L,  and  C with  G - 0. 

These  values  are  R = Rv,  L = Ly,  and  C = = yCa,  for  any  line  of 

constant  axial  cross  section.  Using  these  values  in  the  defining 
equations  of  the  series  impedance  and  shunt  admittance  [equations  (28) 
and  (29)],  together  with  equations  (32)  and  (34),  the  low  frequency 
characteristic  impedance  and  propagation  operator  are  easily  determined: 


x = {ip  ■ fw a+5) = (i^>  (88> 


Tl  = ZlyL  ' = «L  + j*L  H 


■yRyCato* 


d+j) 


(89) 
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Therefore,  the  low  frequency  attenuation,  aL,  and  phase  angle  increase. 


6L,  are 


aL  = iT  • ~ • 


6L  = aL 


At  high  frequencies,  the  inductance  and  capacitance  approach  their 
adiabatic  values.  However,  the  resistance  and  conductance  are  frequency 
dependent  and  become  large  with  frequency.  Karam  and  Franke  [9]  have 
developed  approximations  for  the  high  frequency  circuit  components , 
which  are  rewritten  in  terms  of  u>c  as  follows: 


RH  = f La  ^KLmu.c' 


LH  = La  = p’VA 


Gh  = Ca  ^ KLuwc ' (? 

CH  = Ca  = A/yp*  (? 

In  the  high  frequency  region  R <<  wL  and  G <<  wC.  Using  these 
inequalities,  the  high  frequency  attenuation  per  unit  length,  aH,  and 
phase  angle  increase  per  unit  length,  6h>  are  approximately  [11]: 

LhGh  + RhCh 

au  = ; — ($ 

H 2v^T 

eH  = w/lhCh  ' (? 


Substituting  equations  (92)  through  (95)  into  equations  (96)  and  (97), 


aH  * £ + /Rvca“  = ^ [l  + /kL<^c 


(98) 


6y  = U>V. LaCa 


(99) 


It  is  important  to  note  that  the  foregoing  development  is  applicable  to 
any  fluid  transmission  line  of  arbitrary  but  constant  axial  cross 
section. 

From  equations  (90)  and  (96),  it  is  seen  that  both  and  aH  are 
of  the  form  a = Kt^aT,  where  K is  constant,  and  such  a line  has  a rather 
simple  time  domain  solution. 


Time  Domain  Equations 

The  characteristic  time  tc  is  defined  as  the  inverse  of  the 
characteristic  frequency  (tc  = 1/ujc),  and  the  time  domain  is  divided 
into  three  regions  corresponding  to  the  previously  defined  frequency 
regions.  Thus,  the  short  time  region  (t  <<  tc)  corresponds  to  the 
high  frequency  region  (w  >>  u>c)  and  the  long  time  region  (t  >>  tc) 
corresponds  to  the  low  frequency  region  (w  <<  wc>. 

Metzger  and  Vabre  [23]  have  shown  that  when  the  high  frequency 
attenuation  is  of  the  form  = KvaT,  the  step  response  to  a step 
input  U(t)  at  £ inches  down  a semi- infinite  line  is 


prH  = erfc 


1 ( aK2*: 


1/2 


2 Uf(t-T) 


U(t-T) 


(100) 
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where  U(t)  = |Ps|u(t),  u(t)  is  a unit  step  input,  T = l/c  a is  the 
adiabatic  delay  time,  and  f is  the  signal  frequency.  For  long  times, 
Karam  [20]  has  shown  that  the  unit  step  response  may  be  written  as 


P 


rL 


= erfc 


1 

2 


0l2£2\  1/2 
irf  t / 


U(t) 


(101) 


Since  for  long  times  t - t-T,  evidently  the  unit  step  responses  for 
long  and  short  times  are  governed  by  the  same  relation.  It  seems 
reasonable  to  assume  that  the  step  response  will  retain  the  same  form 
in  the  mid- frequency  range.  Thus,  in  general 


It  is  noted  that  this  assumption  avoids  the  transformation  of  the  mid- 
frequency region  into  the  time  domain.  Brown  and  Nelson  [22]  achieved 
some  success  in  the  mid- frequency  range  for  liquid- filled  circular 
lines,  but  the  method  was  numerical  and  was  reportedly  quite  complex  and 
sensitive.  The  assumption  that  the  attenuation  is  of  the  form  a = 
in  the  mid- frequency  range  is  a convenience  which  allows  the  entire 
time  domain  solution  to  be  expressed  as  a single  complementary  error 
function.  The  "constant"  K is  actually  frequency- dependent  in  the  mid- 
range. Karam  [20]  used  the  time  transform  of  the  known  low  and  high 
frequency  asymptotes  obtained  from  the  equivalent  circuit  approximations 
and  connected  these  with  a logarithmic  straight  line  of  the  form 
log10(a/aH)  = A + B logi0  T',  where  T’  = tv/ro2-  This  straight  line  patch 

is  therefore  made  in  the  time  domain  and  any  correspondence  with  the 
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frequency  domain  is  lost  in  the  mid- frequency  range.  Karam  approximated 
the  transient  behavior  of  a/aH  for  circular  lines  as  [20]: 


aH 


1 

for 

t'<  0.0125 

K(t')® 

for 

0.0125  < t < 1. 25  I 

i (103) 

/S7 

for 

t’>  1.25 

1 

where 


m - ^ log 10 


/8?  1 

1-^1 


and  K - (0.0i25) 


-tn 


A qualitative  plot  of  a/aH  versus  t is  given  in  Figure  28,  along  with  a 
representation  of  Nichols'  transformed  solution.  It  is  emphasized  that 
the  dotted  line  of  Figure  28  representing  Nichols'  solution  is  only 
qualitative;  if  the  exact  transform  of  Nichols'  solution  were  known, 
there  would  be  no  need  for  the  straight  line  approximation.  Karam 
defines  tj  = Tv/r02  and  expresses  the  unit  step  response  as 

? f1  * fej  T° 

I 


F-n  = erf{ 


U(t’-t’o) 


(104) 


■his  solution  gives  one  of  the  simplest  representations  of  the  step 
response  of  circuit  lines  and  shoos  very  good  agreement  with  the  much  more 
complicated  analvtioal  and  numerical  solution  of  Broun  and  Kelson  [213 
for  liquid  lines.  Karam' s solution  is  also  in  good  agreement  with  his 
own  experimental  data  and  that  of  Kantola  [2*3  for  air-filled  lines. 

The  pressure  step  response  of  liquid-filled  lines  obtained  using  the 
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Approximation  of  Attenuation  Ratio  for  Circular  Lines 


results  of  Brown  and  Nelson  [22]  and  of  Karaxn  [20]  are  plotted  in 
Figure  29  for  selected  values  of  to-  Karan's  results  are  applicable 
to  liquids  when  y is  unity  and  the  fluid  bulk  modulus  Kg  is  used  to 
compute  the  celerity.  In  [22,24],  th^  nondimens ional  time  is  also 
defined  as  t'=  tv/rQ2.  It  is  seen  from  Figure  29  that  the  predicted 
responses  are  in  good  agreement  over  a wide  range  of  iq. 


Pressure  ( Flow)  Response  of  Noncircular  Lines  to  & 

Step  Input  of  Pressure  (Flow) 

Whereas  Karam  and  others  used  the  quantity  v/rQ2  to  nondimension- 
alise  time,  the  characteristic  time,  tc  = l/wc,  will  be  used  herein. 
Thus , t'  and  t'q  are  redefined  and  a characteristic  time  equal  to  unity 
is  introduced  as  follows: 

t = wct  = t/tc  = &t' 


TO  = wcT  = T/tc  = 6t  q 


1 - W^tc  ~ t^/tc 


(105) 


In  terms  of  the  redefined  nondimens ional  times,  equation  (104)  becomes 

i [x  * *r]  fe]  T° 


Pr  = erfc 


^(t-Xq) 


U(t-to  ) 


(106) 


and  the  expression  for  the  long  and  short  time  step  responses  become 

1/2. 


PrL  « erfc 


i(B;) 


(Kt-tq) 


(107) 
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rison  of  Pressure  (Flow)  Response  to  a Unit  Step  Input  of  Pressure  (Flow) 
iquid-Filled  Circular  Lines 


TH 


1 

C 2 

6( T- T o ) 

» - 

1/2 


U(t-tq) 


(108) 


The  logarithmic  straight  line  approximation  may  be  rewritten  as  follows: 


1 


for  t < 0 . 1 


a 

“H 


KTm  for  0.1  < t < 10 

— for  x > io 

1 + ^ 

0 


(109) 


where 


and 


K = (0.1)"m 


Equations  (106)  and  (109)  approximately  describe  the  step  response 
of  lines  of  constant  axial  cross  section,  either  circular  or  noncircular. 
The  key  point  is  the  definition  and  use  of  wc,  since  the  ratio  w/u'c 
(hence,  t = t/tc)  characterizes  the  response  of  the  selected  transmission 
line.  For  long  times,  say  one  decade  or  more  above  tc  (t  2 10),  the 
line  response  is  characterized  by  predominantly  viscous,  isothermal 
behavior.  At  one  decade  below  tc  (t  < 0.1)  the  line  performance  is 
nearly  inertial  and  adiabatic.  Furthermore,  the  step  response  of  all 
lines  of  the  same  nondimens ional  delay  time  to  = u>cT  may  be  represented 
by  a single  curve  when  plotted  against  t = wct . Equation  (106)  versus 
t is  plotted  in  Figure  30  for  several  values  of  to*  The  fluid  has  been 
assumed  to  be  air  at  80  F (y  = 1.^017).  Figure  30  represents  either 
the  pressure  response  downstream  to  an  upstream  pressure  step  input  or 
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Pressure  (Flow)  Response  to  a Step  Input  of  Pressure  (Flow)  for  Air-Filled 


the  flow  response  downstream  to  an  upstream  flow  step  input,  i.e.,  the 
equivalent  transfer  function  in  the  Laplace  domain  is 


Pr(s)  _ Qr(s)  _ 1 -tr(s) 
Ps(s)  ' Qs(s)  s 6 


(110) 


Pressure  (Flow)  Response  of  Noncircular  Lines  to  a 
Step  Input  of  Flow  (Pressure) 

Two  other  transfer  functions  are  of  interest  in  determining  the 
response  of  semi-infinite  lines.  These  are:  (1)  the  pressure  response 
to  a step  input  of  fiow,  and  (2)  the  flow  response  to  a step  input  of 
pressure.  The  respective  Laplace  transfer  functions  are: 


Pr(s)  _ p 
Q^CsT  “ s 

Qr(s)  _ 2 

P j (si  s 


L0(s ) 

Y0(s) 


-*r(s) 

e 

-U(s) 

e 


(111) 


(112) 


and  it  is  noted  that  a convolution  integral  may  be  used  to  solve  for  the 
time  domain  response.  Karam  [20]  has  used  the  following  equations  to 
approximate  these  responses : 


■T  p p 

Z0U)  (t-A)dA  ? ^ (t) 
0 s s 

Y0(A)  ^ (t-A)dA  S (t) 
n fs  fs 


r. 


i: 


Z0( A )dA 


Y0(A)dA 


(113) 


(ll1*) 
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The  approximations  given  by  the  terms  on  the  right  hand  side  of 
equations  (113)  and  (114)  are  made  to  simplify  the  results  by  avoiding 
computations  of  the  convolution  integral.  These  approximations  are 
based  on  the  observation  that  Pr/Ps(t)  rises  fairly  rapidly  after  the 
nominal  time  delay  T and  is  essentially  constant  thereafter,  especially 
for  small  t,  whereas  Z0(YQ)  is  essentially  constant  for  small  t,  and 
varies  inversely  (directly)  as  /t  for  large  t.  It  can  be  seen  in 
Tigure  30  that  for  large  tp,  Pr/Ps(t)  does  not  rise  rapidly.  Thus  the 
validity  of  the  approximations  decreases  with  increasing  to*  Fortunately, 
lines  of  large  tp  are  not  practical  transmission  lines;  this  conclusion 
will  be  substantiated  and  illustrated  in  a subsequent  discussion. 

Using  the  large  and  small  s approximations  of  ZD(s)  and  Y0(s),  corresponding 
to  short  and  long  time,  respectively. 


The  integrals  on  the  right  hand  side  of  equations  (113)  and  (114)  mav 
be  easily  evaluated  using  the  inverse  Laplace  transform  [L  ^3  and  the 
approximations  of  equations  (.115)  and  (lib).  Thus 


(life) 


~oH^>dX 


■■so 


•t 

' 0 


\OL(X)d\  = L~l 


1 SO 


(119) 


' 'oH(x)dx  = r‘  [)f  ] ■ X 


SO 


(120) 


For  t ' = tv /r0‘  and  t = w^t  equations  (11)  and  (119)  become 


-oL('V>dX  = “SO 


(W  ----- 


(121) 


~ ^ so  ^ 


S^V  = Yso 


\ ttKlt 


(122) 


It  is  convenient  and  more  realistic  to  construct  a single  continuous 
function  to  describe  the  behavior  of  the  integrals  of  and  10. 
haram  [2C]  has  found  that  equations  (118)  and  (121)  may  be  combined  as 
follows  to  give  a closer  approximation  of  the  overall  behavior  of  the 
integral  of  2C. 


1 


t 

0 


Z0U)dX 


(123) 


Similarly  a continuous  function  approximating  the  integral  of  is 
obtained  in  terms  of  x'  and  t as 
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p 


(124) 


Equations  (123)  and  (124)  versus  t'  are  plotted  in  Figure  31,  along  with 
the  asymptotic  limits  given  by  equations  (116),  (120),  (121),  and  (122), 
for  air- filled  circular  lines  with  = 1.4017.  It  is  seen  that  the 
two  continuous  functions  of  equations  (123)  and  (124)  are  properly 
bounded  by  the  respective  asymptotes.  Using  equations  (113),  (114), 
(123),  and  (124)  the  desired  step  responses  are  obtained  in  dimension- 
less form,  as 


.1/2  i/: 

(t’)  ( 1 + =^(  t)  (l  + — ^ 

"so  \V!s  / fs  \ / rs  V / 


(125) 


'Or 

"so  1 p 
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— (t'  ) 
Fs 


1 t 


8tit 


-1/2 


(t)  1 + 


Kl^t 


-1/2 


(126) 


Equations  (125)  and  (126)  are  plotted  in  Figures  32  and  33,  respectively, 
along  with  the  results  of  Brown  and  Nelson  [22]  for  liquid-filled  lines 
(y  = 1)  and  t'  = tv/r0*-.  It  is  seen  that  the  two  approaches  are  in 
good  agreement  over  a practical  range  of  The  above  step  responses 

are  adapted  to  non-circular  air-filled  lines  as  was  done  for  Fr/Fs  , and 
the  results  are  shown  in  Figures  34  and  35.  In  these  latter  figures, 
the  dimensionless  times  are  based  on  the  characteristic  frequency, 
i.e. , t = wct. 
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Section 


Physical  Implications  of_  the  Fluid  Line  Step 
Responses  and  Dimensionless  Parameters 

It  is  of  interest  to  briefly  discuss  the  physical  implications  of 

the  various  step  responses  and  dimensionless  parameters  presented  in 

this  section.  The  responses  apply  to  semi-infinite  lines,  and 

therefore  represent  waves  travelling  in  one  direction  only,  with  no 

reflections.  However,  the  basic  linearity  of  the  results  allows 

superposition  and  application  to  problems  involving  reflections,  as 

illustrated  by  Karam  [20]  and  Brown  and  Nelson  [22].  The  flow  response 

to  a pressure  step  input  must  eventually  reach  zero  (no  flow)  since  the 

impedance  of  the  semi-infinite  line  is  infinite.  Similarly  the 

pressure  response  to  a flow  step  input  must  eventually  increase  until 

the  line  is  ruptured.  For  lines  having  the  same  tq,  the  responses 

will  be  the  same  when  plotted  against  t.  However,  the  actual  delay 

time,  T,  and  rise  time  of  the  response  will  be  shorter  for  a line 

cross  section  of  higher  u>c.  In  other  words,  lines  of  the  same  1q 

working  fluid  properties  but  different  wc  necessarily  imply  different 

axial  locations  (lengths)  along  the  semi- inf inite  line.  It  follows 

that  the  distance  along  the  line  will  be  shorter  for  the  line  of  higher 

u'c.  The  line  of  higher  uic  will  also  attenuate  high  frequencies  at  a 

more  rapid  rate,  and  will  approach  the  behavior  of  the  diffusive  RC 

line  more  rapidly.  For  signal  transmission,  lines  of  higher  wc  should 

be  avoided,  as  these  are  accompanied  by  larger  attenuation  per  line 

wavelength  at  any  physical  frequency. 

As  to  becomes  large,  say  to  2 1,  the  line  becomes  essentially  a 

diffusive  RC  line,  with  R = Rv  and  C = Cj.  For  a circular  line  of 


1/8  inch  inside  radius  filled  with  air  at  80°F  (Table  I properties) 
and  To  = wcT  = 1,  the  corresponding  length  is  122  ft  and  T is  0.107  sec. 
For  the  same  line  cross  section  and  tq,  but  filled  with  water  at  80°F, 
the  corresponding  length  is  9,576  feet  and  T is  1.945  seconds.  It  is 
unlikely  that  signal  transmission  would  be  of  interest  over  such  long 
fluid  lines.  For  such  lines,  treatment  as  a lumped  RC  line  would  be 
sufficient  for  most  purposes. 

Summary  of  Transient  Responses 

In  this  section  a method  was  presented  for  determining  the 
transient  response  of  semi-infinite  fluid-filled  linesof  arbitrary  but 
constant  axial  cross  section.  The  results  were  presented  in  terms  of  a 
characteristic  frequency  which  depends  only  on  the  laminar  steady  flow 
resistance  and  conductance  of  the  arbitrary  line.  Equivalent  circuit 
parameters  (R,  G,  L and  C)  were  developed  and  used  to  obtain  the  time 
transforms  of  the  various  combinations  of  pressure  and  flow  responses 
to  step  inputs  of  pressure  and  flow.  Results  were  compared  with  those 
of  others  for  both  compressible  and  incompressible  fluids.  For  non- 
circular  lines  the  results  are  shown  for  air  at  80°F.  It  is  noted 
that  the  results  are  restricted  to  laminar  flow.  Similar  results  for 
the  various  step  inputs  propagating  through  fully  developed  turbulent 
mean  flow  are  presented  in  Section  VI. 


93 


..  . . ■i-Ts, 


V . Laminar  Oscillatory  Flow  Combined  with  Steady 
Developed  Mean  Flow 


Background 

Several  investigations  have  been  completed  concerning  the  trans- 
mission of  small  oscillatory  signals  through  fully  developed  mean  flow. 

A majority  of  these  have  been  directed  toward  the  superposition  of 
laminar  oscillatory  signals  on  fully  developed  laminar  mean  flow  [2,3,4]. 
Recently,  investigations  concerning  the  behavior  of  oscillatory  signals 
combined  with  fully  developed  turbulent  flow  have  been  completed  [15, 
25,27].  Most  of  the  recent  efforts  have  been  restricted  to  incompressible 
adiabatic  flow  in  circular  lines . 

Nichols  [3]  obtained  the  small-signal  frequency  response  of  fluid- 
filled  circular  lines  including  frequency  dependent  dispersion  of  the 
propagation  speed  and  distortion  of  the  oscillatory  signal  wave  shape. 

In  addition,  heat  transfer  effects  were  included.  In  essence,  Nichols 
separated  the  steady  and  oscillatory  fluid  parameters  (as  was  done  in 
Appendix  C)  and  subsequently  solved  the  oscillatory  flow  problem; 
nevertheless,  his  solution  is  applicable  to  flows  combining  steady 
developed  laminar  mean  flow  (Hagen-Poiseuille  flow)  and  laminar 
oscillatory  flow.  Holmboe  and  Rouleau  [28]  have  shown  that  if  the  initial 
conditions  are  assumed  to  be  Hagen-Poiseuille  flow  then  the  governing 
differential  equations  may  be  placed  in  a form  which  is  identical  to 
Nichols’  equations.  The  results  of  Brown  [2],  expressed  in  both  the 
Laplace  and  time  domains,  are  identical  to  those  of  Nichols  for  the 
special  case  of  sinusoidal  excitation  signals. 
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More  recently  Brown,  Kargolis  and  Shaw  [15]  have  obtained  the 
frequency  response  of  liquid-filled  circular  lines  carrying  fully 
developed  turbulent  mean  flow.  The  authors  developed  two  and  three 
region  viscosity  profiles  which  resulted  in  fairly  complicated  solutions. 
It  is  noted  that  at  high  Reynolds  numbers  and  dimensionless  frequencies 
thei'e  were  significant  gaps  in  the  frequency  domain  solutions  wherein 
the  attenuation  of  the  oscillatory  signals  was  not  predicted  by  any  of 
the  authors'  models.  Recently  Funk  and  Wood  [25]  developed  a single 
boundary  layer  model  to  account  for  the  transient  viscous  effects  and 
applied  this  model  to  small  amplitude  sinusoidally  disturbed  turbulent 
flow  in  circular  lines.  The  authors  developed  a transfer  function 
relating  the  local  boundary  layer  velocity  gradient  to  the  core  velocity 
of  the  flow.  The  model  is  comparatively  much  simpler  than  that  of  Brown, 
et  al.  [15];  nevertheless  it  gives  values  of  frequency  dependent  signal 
attenuation  which  are  in  excellent  agreement  with  the  theory  and  the 
experimental  data  of  Reference  [15].  It  is  noted  that  the  utility  of  the 
model  given  in  Reference  [25]  is  limited  to  the  high  frequency  regime  and 
constant  phase  velocity;  Funk  and  Wood  suggest  the  use  of  a quasi-steady 
model  (i.e.  , a constant  LRC  model)  at  low  frequencies.  It  will  be 
subsequently  shown  that  the  model  of  Funk  and  Wood  is  closely  related  to 
the  high  frequency  approximations  developed  by  Nichols  [3]  and  Karam  [20] 
for  laminar  oscillatory  flow  combined  with  laminar  steady  flow,  and  to  the 
concept  of  a characteristic  frequency  and  radius  previously  developed 
herein. 

Trikha  [29]  has  suggested  that  the  frequency  dependent  part  of  the 
friction  in  turbulent  flow  may  be  approximated  as  the  frequency  dependent 


part  of  the  friction  in  laminar  flow  for  a limited  range  of  Reynolds 
number  (in  this  matter,  see  also  [7]). 

In  this  section  a simplified  model  of  combined  steady  turbulent  and 
small  signal  oscillatory  flow  is  developed  which  is  based  on  analysis  of 
the  components  R,  G,  L,  and  C of  the  equivalent  electrical  circuit. 

Figure  2,  and  also  on  the  concept  of  a characteristic  frequency.  A 
method  of  determining  the  frequency  response  of  lines  of  either  circular 
or  noncircular  cross  section  carrying  fully  developed  turbulent  mean  flow 
is  developed.  The  results  are  presented  for  a dimensionless  attenuation 
and  frequency.  For  the  circular  line  model,  results  are  compared  with 
those  given  in  [15,25]. 

In  Section  VI,  Karam's  approximate  method  [19]  is  adapted  to  lines 
of  both  circular  and  noncircular  cross  section  carrying  turbulent  mean 
flow  and  the  pressure  and  flow  transient  responses  of  the  lines  to  step 
inputs  are  estimated.  Before  proceeding  with  the  broader  purposes  of 
this  section,  it  is  instructive  to  review  the  physical  aspects  of 
combined  oscillatory  flow  and  steady,  developed  mean  flow  of  compressible 
or  incompressible  fluids  in  a semi-infinite  circular  line. 

Physical  Considerations  of  Combined  Oscillatory  and 
Steady  Developed  Mean  Flow 

A large  body  of  literature  exists  concerning  steady,  fully  developed, 
laminar  or  turbulent  flow  in  circular  lines,  from  which  one  may  select 
flow  relationships  ranging  in  complexity  from  relatively  straightforward 
to  very  complex.  For  turbulent  mean  flow,  elegant  theoretical  flow 
descriptions  involving  statistical  approaches  and  correlation  techniques 
have  been  proposed.  However,  these  descriptions  have  failed  to  yield 


solutions  except  In  the  simplest  cases.  The  approach  taken  herein  is 
to  develop  the  necessary  steady  flow  parameters  from  the  more  useful 
phenomenological  descriptions  of  steady  turbulent  mean  flow.  In  this 
approach,  the  texts  ot  Schlichting  l30j  and  Hint*  131.1  are  particularly 
useful . 

Tor  small  signal  excitation  of  developed  turbulent  flow,  the  steady 
flow  velocity  profile  will  not  he  greatly  altered  by  the  superimposed 
oscillatory  signal  provided  the  steady  flow  Reynolds  number  is  suffi- 
ciently large , and  the  oscillatory  signal  amplitude  is  relatively  small. 

Steady  and  Osci 1 latory  boundary  1. avers 

Tor  steady  developed  laminar  flow  in  a circular  line  or  duct  the 
boundary  layer  is  presumed  to  extend  from  the  wall  to  the  centerline  ot 
the  circular  cross  section.  Thus  the  steady  flow  laminar  boundary  layer 
thickness  denoted  , is  always  equal  to  the  radius  of  the  circular 
line  fiSV{  b r0 1 , and  is  not  dependent  on  the  laminar  flow  Reynolds 
number.  The  classical  Hagen-l'ovseni  lie  velocity  profile  is  obtained 
tor  all  laminar  flow  in  circular  lines.  However , tor  steady  developed 
turbulent  flow  the  associated  velocity  profile  and  boundary  .layer  thick 
ness  will  depend  on  the  Reynolds  number.  As  the  turbulent  ! low  Reynolds', 
number  increases  , the  associated  velocity  profile  changes,  shape  and 
approaches  a plug-flow  profile.  Simultaneously  the  associated  turbulent 
boundary  layer  thickness  Ay ^ decreases,  indicating  large  velocity 
gradient:  at  or  near  the  wall  and  large  resistive  losses.. 

for  laminar  oscillatory  flow,  the  dynamic  boundary  layer  thickness 
will  vary  with  frequency.  Nichols  1 3 1 proposed  a boundary  layer  thick- 
ness, t\„  to,  for  circular  lines;  and  it  is  seen  that  the  oscillatory 


boundary  layer  thickness  decreases  with  increasing  frequency.  Similar 
effects  occur  in  electrical  conductors  where  is  called  the  skin  depth. 
At  low  frequencies  the  oscillators'  velocity  profile  approaches  the 
classical  Hagen-Poiseuille  distribution. 

When  laminar  oscillatory  flow  is  combined  with  steady  mean  flow, 
it  is  assumed  that  an  interaction  between  the  steady  flow  boundary  layer 
and  the  dynamic  (oscillatory)  boundary  layer  will  occur.  If  the  steady 
flow  boundary  layer  thickness  is  less  than  that  of  the  oscillatory  flow 
(i.e.,  <5vt  < or  5Vt  < 6W)  then  the  oscillatory  flow  will  be  dependent 
upon  the  steady  flow  profile.  On  the  other  hand,  when  Sw  < 6Vj,  (or 
«w  < 6yt)  the  behavior  of  the  oscillatory  flow  tends  to  become  independent 
of  the  mean  flow.  This  observation  is  substantiated  both  analytically 
and  experimentally  by  the  work  of  Brown,  et_  al.  [15],  Franke,  karam  and 
Lvmburner  [14],  and  others  [9,25].  Stated  in  other  terms,  the  behavior 
of  the  oscillatory  signal  in  the  low  frequency  region  depends  on  the 
mean  flow  whereas  in  the  high  frequency  region  the  behavior  of  the 
oscillatory  signal  becomes  essentially  independent  of  the  mean  flow  and 
dependent  on  the  oscillatory  signal  frequency  (w).  For  the  experimental 
data  available,  the  high  frequency  behavior  appears  to  approach  laminar 
oscillatory  flow  behavior  [15,27].  It  is  further  assumed  that  for  any 
steady  flow  there  exists  a frequency  at  which  the  dynamic  and  steady  flow 
boundary  layer  thicknesses  are  comparable  ( 6U,  « 6V).  In  this  region,  a 
transition  from  essentially  steady  flow  dependence  to  frequency 
dependence  takes  place,  or  vice  versa.  Similar  qualitative  behavior 
would  be  expected  concerning  the  temperature  profiles,  particularly  for 
fluids  with  a Prandtl  number  near  unity.  For  circular  lines  Nichols  [3] 


defines  a thermal  nonadiabatic  skin  depth  as  67  = /2vj/u  which  is  seen 


to  be  frequency  dependent.  At  low  frequencies  the  heat  transfer  behavior 
of  the  oscillatory  flow  will  depend  essentially  on  the  mean  flow, 
whereas  at  high  frequencies  the  heat  transfer  will  depend  essentially 
upon  the  oscillatory  frequency.  Much  of  the  discussion  of  the  frequency 
domain  equations  in  Section  IV  applies  qualitatively  to  the  above 
discussion.  It  is  not  intended  to  extend  the  discussion  of  the  various 
boundary  layers  any  further,  since  ultimately  a model  for  circular  and 
noncircular  lines  with  mean  turbulent  flow  is  sought.  For  certain  non- 
circular  lines,  more  than  one  5U,  oi'  6V  is  involved  and  it  is  more 
convenient  to  eliminate  the  dependence  of  the  flow  on  the  geometric 
cross  section  of  the  line.  By  integrating  the  various  parameters  over 
the  appropriate  cross  section,  geometry'  dependence  may  be  essentially 
removed,  as  was  previously  done  in  developing  the  characteristic  frequency 
and  components  of  the  equivalent  electrical  circuit  for  laminar 
oscillatory  flow. 

A Mean  Flow  Characteristic  Frequency  for  Circular  Lines 

The  concept  of  a chai’acteristic  frequency  which  divides  the  flow 
into  regions  of  viscous,  isothermal  behavior,  and  inertial,  adiabatic 
behavior  has  been  discussed  in  Section  III  for  circular  and  noncircular 
lines.  It  is  desii-ed  to  extend  this  concept  to  lines  of  arbitrary 
cross  section,  carrying  combined  oscillatory  and  steadv  turbulent  mean 
flow.  Returning  to  the  equivalent  electrical  circuit,  Tigure  2,  and 
keeping  in  mind  that  the  underlying  steady  flow  determines  the  low 
frequency  behavior  of  the  oscillatory  signal,  it  is  necessary  to  determine 
R,  L,  <3,  and  C for  the  underlying  steady  flow.  The  low  frequency 
resistance,  Rv,  is  dependent  on  the  mean  flow. 


(127) 


-3p  -3p 

„ 3z  3z 

Rv  = ~ = ~Jz~ 


and  similarly  for  the  inertance. 


Lv  = KL  ‘ La  = 


i (!) 


dA  • L = 


(128) 


If  the  underlying  mean  flow  is  fully  developed  laminar,  it  was  previously 


shown  that  for  circular  lines 


8tqj 


Lvt  = 3 La 


6v 

“<* = zj 


where  subscript  £ denotes  laminar  flow  conditions.  These  parameters  do 
not  depend  on  Reynolds  number.  For  turbulent  mean  flow,  the  shape  of  the 
velocity  profile  is  a function  of  the  friction  factor  and  Reynolds 
number.  The  defining  relation  for  the  steady  flow  friction  factor,  f, 
for  laminar  or  turbulent  mean  flow  in  circular  or  noncircular  lines,  is 


-3p  f 1 -•> 

~3l  = 2 PU~ 


(129) 


From  equations  (127)  and  (129),  the  laminar  or  turbulent  steady  flow 
resistance  of  any  line  of  arbitrary  but  constant  cross  section  is 


fpu  _ (fRe)u 


(130) 


The  simplest  form  of  the  relation  between  f and  Re  for  developed  turbulent 
flow  in  both  circular  and  noncircular  lines  is  the  Blasius  relation. 


f = a 3164/Re 


(131) 


I 


h 


where  the  Reynolds  number.  Re  is  defined  in  terms  of  the  hydraulic 
diameter: 

Re  = SSSh  (132) 

Equation  (131)  gives  good  accuracy  up  to  Re  < 5*105  [30],  and  is  used  to 
compute  f herein.  Using  Equations  (130)  and  (131)  the  turbulent  viscous 
resistance  of  lines  of  arbitrary  cross  section,  KVt , becomes 


_ 0j_316j4y_  (Re)3/*. 
1 2ADh‘ 


(133) 


From  Hinze  [31],  the  steady  turbulent  velocity  profile  in  circular  lines 
may  be  expressed  as  a simple  power  law. 


* = 

umax  \*o/ 


1/n 


(13«0 


where  n is  a function  of  Re , y = rQ  - r,  and  umax  is  the  maximum 
Icenterline)  velocity.  Further,  the  work  of  Runner  as  presented  bv 
Hinze  [31]  shows  that  the  exponent  n may  be  expressed  as  a function  of 
f only , 


±.  St 

n 


(135) 


Equation  (135)  is  shown  to  be  in  good  agreement  with  the  experimental 
data  of  Nunnor,  Nikuradse,  and  Laufer,  for  both  smooth  and  rough  circular 
lines  [31].  From  Schlichting  [30], 


3n* 


umax 


( n+ 1 ) ( 2n+ 1 ) 


(136) 
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Thus 


u 

u 


(n+l)(2n+l) 


(137) 


The  turbulent  viscous  inertance,  Lvt , becomes 


(n+l)(2n+l)2 

4n2(n+2) 


La  = LLt  • La  (136) 


and  KLt  is  seen  to  be  a simple  expression.  In  terms  of  the  friction 
factor 


KLt 


Lvt  _ {Sf  t l)(/f  t 2)2 
La  U(2^f  + 1) 


(139) 


The  value  of  Klt  will  be  near  unity  since  the  steady  turbulent  velocity 
profiles  are  similar  to  inertial  or  plug  flow  profiles.  As  Re  becomes 
large,  it  can  be  seen  that  u/u  and  kLt  approach  unity,  while  n 00  and 
f -*  0.  It  is  further  noted  that  KLt  near  unity  implies  that  Lv-t  does 
not  differ  much  from  La,  and  the  assumption  that  the  frequency  dependent 
inertia  is  in  fact  constant  (L  = La,  or  L = LVt ) is  a good  approximation, 
particularly  at  high  Reynolds  numbers.  The  approximation  LVt  = La  is 
particularly  useful  in  dealing  with  lines  of  arbitrary  cross  section 
carrying  turbulent  through  flow,  where  expressions  for  u/u  are  either 
very  complex  or  unknown. 

A broad  range  of  values  of  for  circular  lines , together  with 
closely  related  parameters,  were  obtained  using  equations  (131),  (.135) 
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and  113b),  and  are  included  in  Table  VI.  It  is  seen  that  the  tabulated 
uata  is  consistent  with  the  preceding  discussion  of  . 


TABLE  VI 

Turi ulent  Steady  flow  Inert ance  Factor  ) for  Circular  Lines 


Re 

f 

fRe 

n 

1,000 

0.056265 

56.3 

4.2158 

1.0500 

5,000 

0.03762"' 

168.1 

5.1553 

1.0352 

10,000 

0.031640 

316 . 

5.6219 

1.0302 

50,000 

0.021159 

1,057.9 

6.674"1 

1.0211 

100,000 

0.017792 

1,779.2 

7.4969 

1.0180 

500,000 

0.011899 

5,949.3 

9.1675 

1.0125 

1,000,000 

0.010005 

10,005.0 

9.9973 

1.0106 

Characteristic  Frequency  and  Radius  - Turbulent  Mean  Flow 

Proceeding  as  in  Section  III,  a set  of  characteristic  frequencies 
and  radii  based  on  turbulent  mean  flow  quantities  ai'e  defined  as 


Kvt  A|  tv" 

^t  = ~ ; r^t  = ^ ^7 

Rvt  A | 'Sv' 

“Vt  = IT"  ' rVt  " 1 -vt 


(140) 

(141) 


where  subscript  t denotes  turbulent  mean  quantities. 

It  is  more  convenient  to  use  the  characteristic  quantities  defined 
in  equation  (141),  since  these  quantities  (i.e.,  adiabatic  inertance) 
are  more  easily  determined,  and  are  consistent  with  the  previous 
discussion  and  other  publications  on  the  subject  [3,25].  Both  wc^ 
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and  are  dependent  on  the  turbulent  friction  factor  and  Reynolds 

number.  In  essence,  the  circular  line  carrying  fully  developed  turbulent 
flow  is  replaced  by  an  equivalent  smaller  line  of  radius  rv  or  rCt 
carrying  fully  developed  laminar  flow,  of  the  same  characteristic  frequency . 

Nichols'  frequency  domain  solutions  [3],  applicable  to  circular 
fluid  lines  carrying  either  developed  laminar  mean  flow  or  no  mean  flow, 
are  to  be  used  to  solve  the  normalized  frequency-dependent  performance 
of  both  circular  and  noncircular  lines  carrying  fully  developed  turbulent 
flow,  and  are  repeated  herein  for  the  convenience  of  the  reader.  The 
impedance,  admittance,  propagation  operator  and  characteristic  impedance/ 


admittance  are 


2 = j(*)La  1 - 


Y = ju)Ca  1 + 


5>V-1>  Jip7r°i3/2) 


(1*42) 


(1*43) 


r = /£?  = a ♦ jB 


(1*4*4) 


z0  = */z7y1  = i/v0 


(1*45) 


To  use  the  above  solutions  as  stated,  the  radius  rQ  is  replaced  bv  one 


of  the  characteristic  radii,  equations  (1*40)  or  (141),  and  for 


turbulent  flow  the  further  assumption  v-p  « v is  made.  This  assumption 
implies  a i'randtl  number  of  unity  which  is  thought  to  he  more  repre- 
sentative of  turbulent  processes  130,311.  The  adiabatic  inertance  ami 
capacitance  (hu,  Crt)  are  computed  from  the  actual  line  properties, 
liquations  (141)  through  (144)  are  the  same  equations  which  weiv  used 
to  closely  approximate  the  normalised  frequency  dependence  of  both 
circular  and  noncircwlar  lines,  carrying  either  laminar  or  no  mean  flow, 
and  the  results  were  shown  in  Figures  18,  19,  ?‘t,  and  ?£>.  These 
normalized  results  will  remain  the  same,  when  plotted  against  the 
chat acterist ic  frequency  ratio,  <o/wv. . 

As  mentioned  in  Section  111,  tor  a given  value  of  the  argument  ot 
the  has  sol  functions,  * m/v’  r (or  equivalently,  w/uq.  or  <o/u>v),  t ht* 
normalized  frequency  dependence  will  t>e  the  same  when  plotted  against 
this  argument.  Since  it  now  is  intended  to  apply  the  above  modeling 
approach  to  circular  and  nonc.iroular  lines  carrying  turbulent  through 
flow,  it  is  necessary  to  closely  examine  the  range  ot  validity  and 
application  ot  this  approach.  The  following  development  will  address 
the  low,  mid,  and  high  frequency  domains  relative  to  the  turbulent 
viscous  characteristic  frequency,  u>x.  , and  Nichols’  solutions. 

how  Frequency  homa i n - Turbulent  Mean  Flow . For  m ««  a 
quasi-steady  constant  t.RC  model  is  the  simplest  approach,  as  suggested 
by  brown,  ot  aj.  tlh]  and  Funk  amt  Wood  (3(0.  The  equivalent  electrical 
circuit  is  comprised  of  the  constant  value  elements,  K * KVj , 

l,^(*  l,v,  ) » F,  - 0,  and  0 = C\  e \C(X,  The  low  frequency  series 
Impedance  and  shunt  admittance  are 


ZLt  = Rvt  + D^a  = RVt  [l  + j = >La  ^ ~ 3 


YLt  = G + D“ci  - jo»YC£ 


(146) 


(147) 


From  equations  (141)  and  (146),  it  is  seen  that  as  w approaches  zero, 
ZLp  approaches  the  steady  turbulent  DC  resistance,  RVt,  equation  (133). 
The  low  frequency  propagation  operator  and  characteristic  impedance  are 


TLt  = •'ZLtYLt'  = «Lt  + jSLt  = ^ [v  (l  ' 3 ““)]  ( 

-=5 &(“»?)]“■ 

For  very  low  frequencies,  a)/<joVt  <<  1,  the  above  expressions  may  be 
further  simplified: 


(148) 


(149) 


rLt  * \jYRvtCaw  = — ^ jYu>vtw 


Kvt  . fa 
jya)Ca  A u) 


= JYX 

ca  \ 2 “vt 


(1-j) 


(1+j) 


(150) 


(151) 


where  p — ■ » Zs  r ^La/Ca. 


To  facilitate  discussion  of  the  constant  LRC  model  and  comparison  with 
recent  work,  the  dimensionless  frequency,  ft,  and  normalized  attenua- 
tion, aBr>  f°r  circular  lines  introduced  by  Brown,  et_  al^  [15]  will  be 
used.  These  are 
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(112) 


abr  = ~7“  ° 


(153) 


Using  the  constant  LRC  model,  equation  (148),  and  the  above  parameters, 
it  can  lie  seen  in  figure  36  that  a principal  effect  of  increasing  the 
steady  turbulent  Reynolds  number  is  to  increase  the  attenuation  at  all 
frequencies,  up  to  the  laminar  reference  curve.  Tigure  3D  shows  the 
extension  of  the  constant  LRC  curves  beyond  the  laminar  limit,  and 
also  indicates  that  for  a given  Reynolds  number  a plateau  is  reached, 
after  which  the  LRC  model  fails  to  predict  further  increases  in 
attenuation. 

To  obtain  the  limiting  value  of  a,  denoted  a^,  for  the  constant 
LRC  model  and  large  values  of  u>,  rpt  is  first  rewritten  using 
equations  (14D)  and  (147),  and  then  expressed  in  terms  of  its  modulus  and 
phase  angle. 


rLt  = N^t^t  = 


— RveC  £ 


1 - ) 


W\>t 


(154) 


rLt  = 4>*LaC 


T [i  * /t  [•  - <-«*"  '¥] 


(155) 


J'oi  large  values  of  u)/a>vt,  or  small  uiy^/to. 


[■  • (?)’] 
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Variation  of  Low  Frequency  Attenuation  with  Reynolds  Number  in  Circular  Lines 


Taking  aM  as  the  real  part  of  TLt.  and  using  the  relations  RVt  = 
C i = '"Ca*  and  ^CaLa  = l/ca» 


(156) 


Using  the  small  angle  approximation,  sin6  “ 6, 


“M 


1 ^t  I . 1 “vt  n 

2 m 2 ca 


(157) 


Similarly,  for  the  plateau  phase  la'g,  B^, 


(156) 


The  asymptotes  of  value  aM  for  each  Reynolds  number  are  also  plotted 
in  figure  36,  and  it  is  seen  that  aM  is  an  accurate  representation  of 
the  constant  LRC  attenuation  in  the  mid-frequency  range.  Equating  the 
expressions  for  o given  by  equations  (150)  and  (15,1,  the  low  frequenev 
corner,  or  break  frequency  ratio,  u.t/wvt  = 1/2,  is  obtained. 

To  convert  to  Brown's  dimensionless  frequency,  used  in  Figure  36, 


( isy ) 


The  locus  of  break  frequenev  ratios  in  terms  of  ft  is 


I 


This  locus  is  also  plotted  in  Figure  36.  An  expression  for  the 
dimensionless  frequency  fi  at  which  the  asymptote  intersects 
Nichols'  laminar  reference  curve.  Figure  36,  is  obtained  by  equating 
and  a^,  equations  (157)  and  (98),  respectively.  Equation  (98)  was 
derived  from  laminar  theory,  and  contains  developed  laminar  flow 
quantities,  denoted  by  subscript  1. 


1 ^t 

2 Ca 


(161) 


which  can  be  rearranged  to  give 


(162) 


For  incompressible  flow,  y = 1,  and  equation  (162)  becomes 


(163) 


In  terms  of  Brown's  dimensionless  frequency,  see  Figure  36, 


1 


The  entire  set  of  constant  LRC  model  attenuation  curves.  Figure  3n, 
can  be  reduced  exactly  to  a single  curve.  Figure  37,  valid  for  all 
Reynolds  numbers  and  line  cross  sections  by  using  the  dimensionless 
frequency  ratio  w/wv^.  and  expressing  the  dimensionless  attenuation  as 
t'v^ca 

«Lt  = Pv  aLt  {16e) 

where  is  the  real  part  of  rpt  in  equation  (150). 

It  is  seen  from  equations  (150),  (lS^),  and  (162),  and  shown  in 
Figure  37,  that  asymptotic  parameters  are  principally 

functions  of  or  Equation  (166),  with  determined  from 

equation  (146)  and  K = v'  y , is  plotted  in  Figure  37,  together  with  the 
appropriate  asymptotes  and  corner  frequencies.  A second  curve,  which 
has  been  normalized  such  that  the  plateau  value  of  dimensionless 

attenuation  is  unity  (k  * 4>  y"  ),  is  also  shown  in  Figure  37.  It  can 

\ 

be  shown  via  considerable  algebra  that  this  curve  is  also  identical 
to  the  ratio  »y  c/ca  for  *y  X/\a)  obtained  when  using  the  constant 
LKC  model . 

Hid- Frequency  Foma in- Turbulent  Mean  Flow 

If  rv  is  used  in  place  of  rc  in  equations  (142)  through  (144), 
Nichols’  circular  line  theory  predicts  the  low  frequency  attenuation 
very  accurately  at  or  below  the  low  frequency  break  (ti>/wv  = 1 2), 
and  can  be  adjusted  to  predict  the  high  frequency  parameters  also. 
However,  in  the  mid-froquencv  range,  u*  ° u\^,  Nichols’  theory  does  not 
predict  the  extent  of  the  attenuation  plateau;  in  this  region  the 
constant  ERC  model  is  used  and  the  plateau  extended  approximately  to 
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Figure  37.  Attenuation  and  Phase  Velocity  in  Circular  or  Noncircular  Lines,  Laminar  or  Turbul 
H»an  Flow,  Constant  LPC  Model 


the  high  frequency  laminar  reference.  Accordingly,  the  results 
plotted  in  figures  36  and  37  and  the  equations  describing  the  constant 
L.RC  model  are  extended  to  include  the  mid-frequency  range.  This  model 
is  easily  understood  and  applied,  and  is  in  sufficient  agreement  with 
experimental  data  [15,27,32].  Furthermore,  in  the  region  of  the  inter- 
section of  the  constant  LKC  model  plateau  and  the  high  frequency 
laminar  reference  (called  the  high  corner  or  high  frequency  break), 
experimental  anomalies  have  been  encountered  which  are  net  amenable 
to  linear  analysis.  In  a recent  paper  by  Margolis  and  Brown  [27], 
unexpected  and  dramatic  variations  in  both  attenuation  and  phase 
velocity  were  encountered  in  the  transition  region  between  the  quasi- 
steady (constant  LRC),  and  high  frequency  (laminar  oscillatory ) regions 
of  flow  in  water-filled  circular  lines  carrying  turbulent  mean  flow. 

This  same  behavior  was  also  observed  by  Vining  [32]  in  a recent 
experimental  investigation  involving  air-filled  rectangular  lines 
carrying  turbulent  mean  flow.  Vining  and  the  author  conducted  a thorough 
search  for  experimental  anomalies  or  errors  to  no  avail.  The  data  was 
repeatable  from  day  to  day,  after  extensive  experimental  equipment 
replacement  and  recalibration.  A similar  exhaustive  seai'ch  was  made  by 
Margolis  and  Brown,  and  they  were  eventually  persuaded  that  the  experi- 
mental data  was  correct.  Details  of  the  cited  behavior  and  discussion 
of  the  possible  causes  are  best  given  in  [27,32].  However,  the 
behavior  is  worthy  of  illustration,  see  figures  38  and  39.  It  is 
strongly  suspected  that  resonance  phenomena  involving  coupling  between 
turbulent  vortices  and  the  superposed  oscillatory  signal  is  present. 

In  [27],  and  as  seen  in  figures  38  and  39,  the  experimental  data 
corroborates  the  use  of  a quasi -steady  model  for  dimensionless  frequency 
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Laminar 


Constant  IRC  Model 
Laminar  Reference  [3] 
Data  of  Margolis 
and  Brown  [27] 


Normalized  Attenuation  Data  in  Water  Medium 


Constant  IRC  Model 
Laminar  Reference  [3] 
Data  of  Margolis 
and  Brown  [27] 


Figure  39.  Phase  Velocity  Ratio  Data  in  Water  Medium 


ratios  below  and  up  to  the  transition  region,  and  the  use  of  a laminar 
oscillatory  model  at  ratios  above  this  region.  The  experimental  evi- 
dence in  f 27, 32]  further  emphasizes  the  futility  of  attempting  to 
develop  a rigorous,  complicated,  smooth  transition  model  as  was  done 
in  [15,25j;  attenuation  and  phase  velocity  behavior  in  the  transition 
region  is  tar  from  smooth,  fortunately,  the  transition  frequency 
bandwidth  is  predictable  [27],  and  can  be  avoided  in  any  attempted 
transmission  of  information  through  the  developed  turbulent  mean  flow. 

In  view  of  the  experimental  evidence,  it  is  thought  sufficient  to  use 
the  constant  LRC  model  for  dimensionless  frequency  ratios  in  the  mid- 
range and  up  to  those  at  which  the  laminar  oscillatory  model  becomes 
valid. 

High  Frequency  Domain  - Turbulent  Mean  flow.  In  the  high  frequency 
domain  (m  >>  mv^),  the  flow  becomes  frequency  dependent.  At  high 
frequency,  the  laminar  oscillatory  boundary  layer  becomes  considerably 
thinner  than  the  equivalent  mean  flow  boundary  layer  (6fa,  <<  ) as 

previously  discussed  in  this  section.  Thus  the  high  frequency  behavici 
becomes  essentially  independent  of  the  mean  flow,  and  approaches  the 
laminar  oscillatory  behavior  associated  with  the  line  being  investigated. 
The  normalized  laminar  oscillatory  behavior  of  both  circular  and 
noncircular  lines  was  shown  to  be  nearly  identical  when  plotted  against 
the  laminar  characteristic  frequency  ratio  u>v^  = Rvf/i-a*  remains  to 
determine  the  transition  from  middle  to  high  frequency  and  to  normalize 
the  high  frequency  domain  results  such  that  they  are  applicab  e to 
circular  and  noncircular  lines  carrying  laminar  or  turbulent  mean  flow 
over  a practical  range  of  Reynolds  number.  Results  are  sought  which 
are  similar  to  the  single  curves  of  Figure  37.  These  curves  include 


both  low  and  mid-frequency  performance  of  circular  and  noncircular 


linos  carrying  a compressible  or  incompressible  fluid  at  the  speciried 


Reynolds  number. 


If  Nichols'  equations  (142  through  145)  are  used  together  with  rVt 
the  predicted  high  frequency  attenuation  is  too  high,  as  shown  in 
Figure  40.  Thus  rv  cannot  be  used  as  the  characteristic  radius  in  the 
high  frequency  region  unless  other  adjustments  are  made.  Experimental 
evidence  supports  the  use  of  the  laminar  oscillatory  model  in  the  high 
frequency  regions,  and  the  following  development  will  provide  the 
mathematical  relationships  and  adjustments  needed  in  the  high  frequency 
domain. 

Equation  (144)  may  be  expanded  using  equations  (142)  and  (143) 
with  v = \>j  to  give 
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The  normalised  impedance  ratio  is  obtained  from  equations  (142), 
(145),  and  (36): 
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It  is  clear  that  the  argument  <f>  completely  determines  the  bracketed 
expression  of  equati'„i  (167)  and  all  of  equation  (170).  The  modulus 
of  $ can  rewritten  in  terms  of  the  characteristic  radius  and 
frequency  as 


where  rv  = rc,  rv^  or  rVt,  and  similarly  for  o>v.  The  argument  $ depends 
only  on  ot/otv,  or  its  equivalents.  However,  from  equation  (107)  it  is 
seen  that  T also  depends  on  ot/ca.  Supposing  two  lines  of  the  same 
characteristic  frequency  ratio,  but  different  wv,  i.e., 
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Consequently,  a (and  6)  may  be  adjusted  between  lines  of  different 
characteristic  frequency  to  obtain  the  proper  values: 


It  follows  that  a single  baseline  curve  can  be  constructed  to  describe 
the  high  frequency  attenuation  of  all  lines,  carrying  either  laminar 
or  turbulent  mean  flow.  One  form  of  this  curve  is  obtained  by 
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normalising  a.  equation  (107),  in  a manner  similar  to  brown  lit]: 
''0ca  u> 
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when  a = a(rv,).  rquation  (17*0  shows  that  the  normalised  attenuation  i-- 
a function  of  u>/tov  only.  Equation  (174)  is  plotted  in  figure  41.  The 
high  frequency  behavior  as  shown  in  figure  41  (above  to/uq^  = 7),  may  be 
adjusted  to  the  appropriate  laminar  reference  (asymptote)  by  using  the 
Known  high  frequency  behavior  and  t he  approach  leading  to  equations  (l*Vl 
and  (173),  or  by  noting  that  the  correct  characteristic  frequency  ratio 
at  high  frequencies  is  to/wVj. . whbther  the  underlying  steady  flow  is 
laminar  or  turbulent,  and  interpreting  figure  4]  accordingly . 

It  is  convenient  to  normalise  u such  that  the  onset  of  frequency 
dependent  behavior  begins  at  unity  value  on  the  mid-frequency  plateau 
and  joins  the  laminar  oscillatory  limit  curve  of  the  fluid  line  being 
considered.  The  normalized  high  frequency  attenuation  is  to  hr  determined 
in  terms  of  u»/u'v^ . The  appropriate  asymptotes  are  given  in  equations 
(1R7)  and  (eg)  for  the  middle  and  high  frequency  regions,  respectively. 
These  are  used  to  define  a new,  extremely  simple,  single  expression  for 
the  middle  and  high  frequency  dimensionless  attenuation.  The  two 
asymptotes  are: 
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Frequency  Dependent  Attenuation  in  Circular  and  Noncircular  Lines,  Laminar  o 
Turbulent  Mean  Flow 
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Equations  (175)  and  (177)  are  combined  into  a single  function,  called 
5Ht: 


“Hr 
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This  function  represents  the  middle  and  high  frequency  dimensionless 
attenuation  of  laminar  oscillatory  signals  in  circular  and  noncircular 
lines,  carrying  fully  developed  laminar  or  turbulent  throughflow.  The 
inclusion  of  y in  this  equation  also  accounts  for  compressibility  effects. 
Equation  (178)  is  plotted  in  Figure  42.  The  high  corner  (break) 
frequency  occurs  at  the  intersection  of  the  two  asymptotes,  as  previously 
shown  in  the  development  of  equations  (161)  through  (165).  A new 
dimensionless  grouping,  the  abscissa  of  Figure  42,  is  formed  which 
combines  mean  flow,  geometry,  and  frequency  dependent  parameters.  This 
new  grouping  facilitates  comparisons  with  previous  work,  notably  that  of 
Funk  and  Wood  [25],  as  discussed  in  the  following  subsection. 

Comparison  of  High  Frequency  Results  with  Previous  Work 

In  comparing  the  above  results  with  those  of  Brown,  et  al.  [15,27], 
and  Funk  and  Wood  [25,26],  it  is  necessary  to  establish  a common  basis 
among  the  references.  This  basis  assumes:  (1)  lines  of  circular 
cross  section  (rV£  = rQ  or  R);  (2)  nearly  incompressible  fluids  (y  = 1); 
(3)  unity  Prandtl  Number  (o  = 1);  (4)  constant  phase  velocity  (c  = ca); 
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(5)  developed  turbulent  mean  flow;  (6)  relatively  small  amplitude 
oscillatory  signals;  (7)  small  attenuation  per  line  wavelength;  and 
(8)  high  or  near  high  signal  frequencies  (w/wVt  > 1),  for  all  turbulent 
Reynolds  numbers.  Since  Funk  and  Wood  [25]  have  compared  their  work  to 
that  of  Brown,  et  al . £153,  the  comparison  herein  deals  primarily  with 
Reference  [25].  The  dimensionless  parameters  of  Brown,  equations  (152) 
and  (153),  will  be  used  for  overall  comparisons,  since  these  parameters 
have  found  wide  use  in  the  literature. 

The  model  developed  by  Funk  and  Wood  [25]  assigns  all  viscous 
effects  to  a single  fictitious  laminar  boundary  layer  whose  thickness 
(A)  is  based  on  steady  flow  laminar  and  turbulent  friction  factor  data: 


A - 16 

R ' (fRe )t 


(179) 


where  R is  the  radius  of  the  circular  line.  All  other  effects  are  assigned 
to  an  inertial  plug-flow  core  (i.e.,  L - La).  The  model  makes  use  of 
the  eight  assumptions  listed  above,  and  results  in  a single  dissipation 
function  ($)  such  that  the  attenuation  can  be  expressed  as 


• « 


' 


v# 

° = RcaA 


(180) 


The  attenuation  itself  is  assumed  to  be  small,  and  only  first  order  eflects 
are  considered  in  the  expansion  of  the  dynamic  pressure  ratio: 


(181) 
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The  model  of  Funk  and  Wood  still  requires  the  computation  of  a complex 
transfer  function  relating  the  local  velocity  gradient  to  the  core 
velocity,  and  also  requires  numerical  integration  of  the  dissipation 
function  across  the  boundary  layer.  Both  of  these  requirements  are 
avoided  by  the  model  given  in  equation  (178)  and  plotted  in  Figure  42. 
From  equations  (65),  (75),  (76),  and  (141)  it  is  easily  shown  that  for 
the  same  circular  line  carrying  turbulent  mean  flow. 
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The  model  of  Funk  and  Wood  and  the  author's  model  are  related  through 
the  turbulent  mean  flow  parameters.  Comparing  equations  (179)  and  (182), 
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The  relationship  between  $ and  the  parameters  developed  herein  is  obtained 
by  comparing  expressions  for  the  dimensionless  attenuation.  From 
equation  (180), 
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and  from  equations  (178)  and  (183), 
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It  follows  that  equivalence  is  of  the  form  & = 5^.  From  equations 
(178),  (183),  and  (185),  with  y = o = 1,  and  rV£  = R for  comparison 
with  [25], 
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Equation  (186)  indicates  the  significant  influence  of  the  steady  and 
dynamic  boundary  layers  on  the  line  response.  At  high  frequency  ratios. 
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and  at  low  frequency  ratios , 
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and  the  high  corner  (break)  frequency  occurs  at  A = 6U,.  From  Nichols  C 3 ] 
and  Karam  [20],  the  high  frequency  laminar  oscillatory  attenuation  was 
approximated  by  equation  (98).  Equation  (98)  reduces  to  equation  (187) 
as  shown  in  the  following  development  (y  = o = 1). 
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Using  the  delinitions  of  u)Vj;  and  6^,  together  with  equation  (189)  and 
Brown's  notation,  gives 
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It  remains  to  show  the  equivalence  of  the  mid-frequency  asymptotes, 
i.e.,  that (R2ca/v)  = R/A,  with  y = o = 1: 


:a  /7  “Vf 
~ * 2 c_ 


(-)' 


(191) 


The  last  two  terms  of  the  above  equation  are  identical  to  those  of 
equation  (183). 

The  attenuation  model  embodied  in  equations  (178)  and  (186), 
applicable  to  the  middle  and  high  frequency  regions,  constitutes  a very 
significant  simplification  when  compared  to  the  models  given  in  [16]  and 
[25].  An  appreciation  of  this  model  can  only  be  realized  by  comparison 
with  the  cited  references. 

The  results  of  Brown,  £t  al.  [15,27],  Funk  and  Wood  [25],  and  the 
author  are  plotted  for  comparison  in  Figure  43  at  two  Reynolds  numbers, 
Re  = lxlO4  and  1x10s,  for  rv^  = R.  In  view  of  the  experimental  data, 
the  author's  model  is  as  valid  as  any  of  the  referenced  models  for 
circular  lines,  and  is  comparatively  much  simpler.  Furthermore,  the 
model  is  applicable  to  any  fluid  line  of  arbitrary  cross  section, 
carrying  laminar  or  turbulent  mean  flow.  The  other  models  are  not. 


Broun,  et  al . , Invariant  Turbulence  Model  [IS] 


bulent  Mean  Flow 


VI . The  Transient  Response  of  Fluid  Lines  Carrying 
Developed  Turbulent  Mean  rlow 


In  this  section,  the  high,  middle,  and  low  frequency  models  developed 
ir.  Section  V are  transformed  into  the  time  domain  using  the  general 
approach  of  haram  [19,20].  The  methodology  and  form  of  the  data  is  very 
similar  to  that  presented  in  Section  IV;  much  of  the  general  discussion 
presented  in  Section  IV  concerning  the  various  regions  of  the  frequency 
and  time  domains  is  also  applicable  here.  Accordingly,  only  major 
differences  will  be  pointed  out  in  the  development  and  discussion  which 
follows.  In  all  developments  of  this  section,  unity  Prandtl  number  is 


assumed . 


Frequency  and  Time  Domain  Equations 


The  frequency  domain  is  divided  into  three  regions  (low,  middle,  and 
high)  based  on  the  turbulent  characteristic  frequency  ratio,  w/u\,t  = 
Rvt/La.  The  regions  are  divided  according  to  the  value  of  the  low  and 
high  corner  (break)  frequency  ratios,  as  defined  by  equations  (150), 


(157),  and  (162): 
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It  is  seen  that  the  regions  are  highly  dependent  on  both  the  laminar 


and  turbulent  steady  flow  parameters.  A new  characteristic  time,  tc  = 
l/wv  , is  used  to  redefine  the  dimensionless  time  parameters,  such  that 

t = u)Vtt  = t/tc 

t0  = wVtT  = T/tc  1 193 ) 

1 — (il  V-.  ^ c *"  tc 

The  above  frequency  regions  and  dimensionless  time  parameters  are 
used  in  developing  expressions  for  the  various  pressure  and  flow 
transient  responses  presented  in  this  section. 

Line  Attenuation  Factor 

To  facilitate  the  use  of  Karam's  methodology,  the  attenuation  in 
turbulent  mean  flow  is  expressed  in  terms  of  the  three  logarithmic  line 
segments  asymptotes ) corresponding  to  the  three  frequency  regions  of 
equations  (192): 
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It  is  seen  that  the  high  and  low  frequency  attenuation  are  of  the  form 
« * h >Tt.  Consequently,  the  solutions  of  Metzger  and  Vabre  [2?]  and 
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Karam  [?0j  may  be  adapted  to  obtain  the  various  transient  responses.  The 
attenuation  ratio,  a/a^,  is  redefined  as 


Pressure  (.Flow)  Response  to  £ Step  Input  of 
Pressure  (.Flow ) , Turbulent  Mean  Flow 

The  general  solution  for  the  transient  downstream  pressure  (flow) 
response  to  a step  input  of  pressure  (flow)  was  given  by  equation  (102). 
Using  the  attenuation  ratios  of  equations  (195)  above,  together  with 
equations  (1C2)  and  (193),  the  generalized  transient  response  may  be 
expressed  as 


(196) 


where  subscript  t denotes  applicability  to  turbulent  mean  flow  conditions. 
The  dimensionless  step  responses  for  long  and  short  times  corresponding 
to  equations  (10?)  and  (10S)  are,  respectively: 
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and 
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For  completeness,  the  intermediate  or  middle  time  solution  is 
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It  is  seen  that  the  middle  and  long  time  transient  response  is  completely 
determined  for  a given  fluid,  T,and  to-  It  is  also  noted  that  equation 
(197)  is  identical  to  equation  (107);  in  the  latter  equation  the  times 
are  nondimensionalized  with  respect  to  i uc.  Equations  (197)  and  (199) 
describe  the  transient  response  of  arbitrary  lines  filled  with  a 
compressible  fluid  carrying  laminar  or  turbulent  mean  flow.  For  laminar 
mean  flow,  wV£  = RVj^/La  is  used  to  define  the  dimensionless  time  (i.e., 
t = u)V£t);  for  turbulent  mean  flow,  u>Vt.  = Rv^/Ea  is  used.  The  short 
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ime  response,  equation  (198),  retains  some  dependence  on  the  ratio  cf 


the  turbulent  to  laminar  characteristic  radii.  The  physical  interpre- 
tation is  that  the  high  corner  break  frequence,  which  essentially  marks 
the  division  of  influence  between  steady  turbulent  mean  flow  behavior 
and  frequency  dependent  laminar  oscillatory  behavior  must  include 
influences  from  both  phenomena.  The  ratio  rv  /rVJ  essentially  adiusts 
the  transient  response  to  account  for  the  various  arbitrary  line 
geometries  and  turbulent  Reynolds  numbers,  and  is  therefore  not  solely 
dependent  on  the  turbulent  characteristic  frequency  (u>yt). 

Equation  (196)  is  plotted  in  Figure  44  as  a function  of  t = u>v^t , 
for  two  values  of  the  ratio  rv, /rV{ , which  correspond  to  the  following 
tabulated  conditions  for  circular  lines  (rv^  = R). 


TABLE  VII 

Equivalent  Circular  Line  Steady  Flow  Parameters  for 
rigures  4_4 , 4b,  and  46 
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.44975 

1.4017 
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316.4 

.15760 

.0935 

. 07998 

1.4017 
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1779.2 

.02803 

.0935 

The  influence  of  rvt/rv^  on  the  short  time  response  is  small  and  nearly 
negligible,  as  can  be  seen  in  Figure  44.  In  the  middle  or  long  time 
regions,  there  is  no  influence,  which  is  evident  from  equations  (197) 
and  (199).  For  the  large  differences  in  mean  flow  parameters  given  in 
Table  VII,  the  nearly  unnoticeable  effect  on  the  pressure  and  flow 
parameters  plotted  in  Figure  44  indicates  that  the  transient  responses 
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are  very  weak  functions  of  rVt/rV£.  The  equations  developed  in  this 
section  also  apply  to  upstream  and  downstream  flow  transients;  i.e., 

K>-(£ 

and 


U(t)  = Ps(t)u(t)  or  Qs ( t )u( t ) 

as  appropriate.  Expressions  for  the  remaining  unit  step  inputs  and 
responses  of  interest  are  developed  in  the  next  subsection. 

Pressure  ( Flow ) Response  to  a Step  Input  of  Flow 
(Pressure),  Turbulent  Mean  Flow 

The  two  remaining  dimensionless  transfer  functions  of  interest  are 


Mf)\  i Zo(s)  ^-ar-t-Cs) 

Qs^s)  / S ‘'so 

't 

and 


(200) 


1 fQr(s)  \ _i  -trt(s) 

T^TsT  y psv s ) J s ^ so 


(201) 


As  was  done  in  Section  IV,  the  approximations  used  bv  karam  [20]  arc 
adapted,  with  r = wv  t,  and  equations  (113)  and  (11*0  are  used  to  obtain 
the  final  time  domain  expressions.  Subscript  t again  denotes  applica- 
bility to  turbulent  mean  flow  conditions. 

Approximations  of  tne  characteristic  impedance  and  admittance  for 
the  long  time  (small  s)  and  short  time  (large  s)  regions,  adapted  to 
turbulent  mean  flow  conditions,  are  respectively : 


13*4 
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:°Ht = w = 2so 


(203) 


Performing  the  required  integrations  yields: 


= I'l 


Zq^  (X)d(X)  = L 
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o 

■t 

V°Ht(X)dX 

= IT1 

>] 

- ^SO 

(207) 

As  before,  it  is  convenient  to  combine  the  impedance  and  admittance  into 
a single  function  applicable  to  the  entire  time  interval: 
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•t 

2c(X)dA  = 
o 


(206) 


•t 

Y0Q)dA  = 
o 


(209) 


Using  the  form  of  equations  (113)  and  (114)  together  with  equations 
(196),  (206),  and  (209),  the  desired  dimensionless  step  responses  are 
obtained  as  follows: 


(210) 


(211) 


Equations  (210)  and  (211)  are  plotted  in  Figures  45  and  46,  respectively, 
for  the  parameters  given  in  Table  VII. 

Physical  Implications  of  the  fluid  Line  Step 
Response  with  Turbulent  Mean  Flow 

It  is  of  interest  to  discuss  the  physical  implications  of  the 

dimensionless  step  responses  shown  in  Figures  44,  45,  and  46.  For  lines 

with  the  same  value  of  tq,  the  computed  responses  will  be  identical 

for  values  of  x i (y/4)  (tuv^/wv^),  and  will  be  nearly  identical  for 
smaller  values  of  t.  However,  two  lines  of  identical  tq  and  working 
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Response  to  a Step  Input  of  Pressure  for  Air-Filled  Lines  of  Arbitrary 
Section,  Carrying  Turbulent  Mean  Flow 


fluid,  but  different  turbulent  characteristic  frequencies  necessarily 
imply  different  physical  lengths  down  the  lines.  The  line  having  the 
larger  turbulent  characteristic  frequency  will  exhibit  a faster  rise 
to  a given  ratio  of  Py/Ps  or  Qr/Qs>  in  a shorter  length,  and  will 
approach  a diffusive  RC  line  more  rapidly.  Comparing  the  two  circular 
pneumatic  lines  for  the  conditions  of  Table  VII  with  tq  = 2,  it  is  seen 
that  an  increase  in  Reynolds  number  from  Re  = 10**  to  Re  - 106  decreases 
the  adiabatic  delay  time  (T  = i/ca)  by  the  ratio  of  31.6/1!  The 
implication  of  this  example  is  that  increases  in  the  turbulent  flow 
Reynolds  number  greatly  accelerate  the  transition  to  diffusive  RC  line 
behavior;  consequently,  high  turbulent  Reynolds  numbers  should  be 
avoided  in  any  attempt  to  transmit  information. 

In  terms  of  the  equivalent  laminar  and  turbulent  steady  flow 
boundary  layers,  equation  (183),  it  is  seen  that  the  condition 


is  equivalent  to 


for  circular  lines  or 

t i 

for  arbitrary  lines.  Thus  the  break  time  between  short  time  and  inter- 
mediate time  domain  responses  is  a strong  function  of  the  equivalent 
steady  flow  boundary  layer  parameters. 


I 
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VII.  Conclusions  and  Recommendations 


This  section  contains  a summary  of  the  conclusions  reached  as  a 
consequence  of  this  work.  Recommendations  for  extensions  of  this  work 
are  addressed. 

Conclusions 

1.  In  general,  the  very  close  analogy  between  the  small-signal 
response  of  fluid  and  electrical  transmission  lines  has  been 
reinforced  for  a wide  range  of  fluid  line  geometries,  mean  flow 
conditions,  fluid  media,  and  heat  transfer  effects.  The  signifi- 
cant role  of  the  equivalent  electrical  transmission  line  properties 
in  describing  the  fluid  line  response  has  been  demonstrated. 
Distributed  line  parameters  such  as  resistance,  inductance, 
conductance  and  capacitance,  derived  from  the  series  impedance 

and  shunt  admittance  of  an  equivalent  electrical  transmission  line, 
are  shown  to  be  the  key  parameters  in  determining  the  line  response. 

2.  A closed  form,  distributed  parameter,  analytical  solution  for  the 
small-signal  frequency  response  of  fluid  transmission  lines  of 
annular  cross  section  has  been  obtained.  The  solution  considers 
shear  flow  effects,  fluid  compressibility,  and  heat  transfer 
effects,  and  is  applicable  to  blocked  lines  or  lines  carrying 
developed  laminar  mean  flow.  The  solution  was  shewn  to  be  in  very 
good  agreement  with  experimental  data  for  blocked  annular  pneumatic 
transmission  lines  over  a wide  range  of  annular  radius  ratios  and 
signal  frequencies.  In  comparing  the  annular  line  results  with 
those  previously  obtained  for  rectangular  lines,  it  was  shown  that 
for  large  radius  and  aspect  ratios  the  solutions  converge.  The 


annular  line  solution  contains  the  circular  line  solution  as  a 


special  case. 

з.  Lines  of  arbitrary  cross  section  can  be  accurately  modeled  as 
equivalent  circular  lines  having  the  same  ratio  of  steady  flow 
resistance  to  inertia  properties.  The  concept  of  an  equivalent 
characteristic  frequency  and  radius  can  be  applied  to  transmission 
lines  of  arbitrary  cross  section  to  predict  both  the  frequency 
response  and  transient  behavior  of  these  lines,  including  either 
laminar  or  turbulent  mean  flow  effects.  In  solving  practical 
problems  involving  fluid  systems,  this  equivalence  should  be  of 
significant  value  in  simplifying  the  system  modeling  task. 

и.  The  characteristic  radius  is  analogous  to  the  hydraulic  radius 
used  to  represent  circular  and  noncircular  lines  in  steady  flow- 
analysis.  However,  the  characteristic  radius  more  accurately 
accounts  for  the  actual  viscous  and  inertial  effects  and  was  shown 
to  be  clearly  superior  to  the  hydraulic  radius  in  modeling  the 
frequency  response  of  noncircular  lines.  In  view  of  the  much  closer 
agreement  with  existing  data  and  analytical  solutions  for  both 
rectangular  and  annular  lines,  the  characteristic  radius  should  be 
used  instead  of  the  hydraulic  radius  to  determine  the  equivalent 
circular  line  for  all  noncircular  lines. 

5.  Very  significant  simplifications  have  been  achieved  regarding  the 
modeling  of  the  frequency  and  temporal  response  of  circular  lines 
carrying  turbulent  mean  flow.  Comparison  of  the  simplified  models 
with  existing  models  and  experimental  data  on  circular  lines  shows 
that  the  simplified  models  are  as  accurate  as  any  existing  models 
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in  predicting  the  frequency  response.  Furthermore,  the  simplified 
models  are  not  restricted  to  circular  lines  or  incompressible 
fluids. 

6.  Of  the  four  equivalent  circuit  components,  the  steady  flow  viscous 
resistance  is  the  most  influential  parameter  in  determining  the 
frequency  and  temporal  response  of  fluid  transmission  lines.  For 
developed  laminar  flow,  the  steady  flow  resistance  is  determined  bv 
the  line  geometry  and  fluid,  and  is  not  influenced  by  changes  in 
the  Reynolds  number.  For  developed  turbulent  mean  flow,  the 
resistance  is  strongly  dependent  on  the  Reynolds  number.  In  all 
cases,  the  thickness  of  the  frequency  dependent  laminar  oscillatory 
boundary  layer  relative  to  the  equivalent  steady  flow  boundary 
layer  is  a key  parameter  in  determining  whether  the  frequency  or 
transient  response  exhibits  predominantly  viscous,  isothermal 
behavior  or  inertial,  adiabatic  behavior. 

Recommendations 

1.  The  small  signal  transient  response  of  both  circular  and  non- 
circular lines  carrying  laminar  throughflow  should  be  determined 
experimentally  and  compared  with  the  analytical  results  presented 
in  this  work.  The  experiment  should  include  both  compressible 

and  incompressible  fluid  flow,  and  should  be  repeated  for  turbulent 
throughflow. 

2.  More  analytical  and  experimental  work  is  also  needed  to  extend  the 
applicability  of  the  concepts.  The  frequency  and  temporal  responses 
of  fluid  transmission  lines  to  large  amplitude  excitation  signals 
are  of  particular  interest.  The  assumption  of  essentially  constant 


mean  fluid  properties  and  small  perturbation  of  a steady  mean  flow, 
valid  for  small-signal  excitation,  would  have  to  be  re-examined. 

A possible  approach  would  begin  with  a sensitivity  analysis  of  the 
amplitude  and  frequency  dependence  in  the  initial  formulation  of 
the  governing  equations , followed  by  the  addition  of  corrective 
terms  to  the  simplified  equations  of  this  work  to  account  for  the 
first  order  dependence  on  excitation  signal  amplitude. 

3.  An  analytical  model  which  adequately  describes  the  behavior  of  the 
signal  attenuation  and  phase  velocity  at  or  near  the  high  frequency 
corner,  or  break  frequency,  should  be  developed.  The  anomalous 
behavior  is  associated  with  the  transition  of  small-signal 
dependence  from  quasi-steady  turbulent  mean  flow  dependence  to 
laminar  oscillatory  frequency  dependence.  A companion  experiment 
which  would  include  flow  visualization  and  high-speed  photography 
of  the  disturbances  might  provide  significant  physical  insights 
in  developing  the  model. 
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Appendix  A 


DC  Resistance  and  DC  Inductance  of  Fluid  Transmission 
Lines  Carrying  Developed  Laminar  Flow 

In  general,  the  DC  resistance  (Rv)  and  DC  inductance  (Lv)  are 
determined  from  the  steady  flow  parameters  of  the  rigid  transmission 
line.  These  parameters  will  be  functions  of  the  line  geometry,  fluid 
properties  and  steady  flow  profiles.  The  DC  resistance  is  defined  by 

_ _ -3P/9z  -3P/3z 

Rv  . _____  . __  U_ 


and  the  DC  inductance  is  defined  as 


(A-2) 


where  u is  the  steady  flow  velocity  (either  laminar  or  turbulent);  u 
is  the  average  velocity,  u = 1/A  udA;  and  La  = p/A  is  the  adiabatic 
inductance. 

For  circular  lines,  Rv  and  Lv  are  constants.  For  steady  developed 
laminar  flow,  Rv  and  Lv  are  obtained  from  the  classical  Hagen-Poiseuille 
flow: 


„ _ -3P/3z  _ Lthj 

Rv  - — 

The  Hagen-Poiseuille  velocity  ratio  is  [30] 


(A- 3) 


(A-4) 
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where  r is  the  radial  coordinate  and  rQ  is  the  inner  radius  of  the 
circular  line.  Using  this  ratio  in  equation  (A-2), 


7TTr 


2irrdr  = ^ La 


(A-5) 


For  annular  and  rectangular  lines  Ry  and  Ly  will  depend  on  the  radius 
or  aspect  ratio,  respectively.  For  the  annular  line  [8], 


re2/4 u (-dP/dz) 


= if1  + r!2  ♦ 


(A-6  ) 


From  which  RVn  is  easily  obtained  as 


- 3P/3z 

Q 


- (*  1+ri2  . 1 1 

" An2  L !-r  • 2 In  r*J 


(A- 7) 


where  r^  = rj/r0,  is  the  radius  ratio  of  the  annular  cross  section. 
The  velocity  ratio  is  [ 8 ]: 


u [1-r*2]  In  rj[  - [1-r*2]  In  r* 
“ ^ C(l+r*2)  In  r*  + (l-r5{2)3 


(A- 8) 


where  r*  = r/r0  and  the  dependence  on  geometry  is  obvious.  The  notation 
r*  = ri/rD  is  generally  used  in  the  body  of  this  work.  However,  for 
this  appendix  only,  it  is  necessary  to  distinguish  between  the 
dimensionless  ratio  r*  = r/r0  and  r^  = ri/r0.  Using  equations  (A-2) 
and  ( A- 8 ) , 
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(A-9) 


where 


d2  = [(i  t ^ ) In  z-t  t (l  - «^)]2 
For  the  rectangular  line  [8]  of  aspect  ratio  a = h/b. 


(A- 10) 


hi 

(-  8P 

4y 

V 82 

= 2 I cos(ain) 

i~l  “i3 


coshl 


cosh(^|) 


(A- 11) 


where  n = y/h,  £.  = x/b,  and  ai  = v/2(2i-l)v/2(2i-l)  i 1,2,3,.. 

The  average  velocity  is 


°°  TaT  ai"l 

2 l ~^7s  sin2°i  ' tanh  ~J 


( A-12 ) 


h2/-3P\  i=l  ai 

4y\  8z  / 

The  DC  resistance  for  the  rectangular  line  is  easily  found  as 

\l  3 

[i=l  ai 


- 8P/8z  _ 2 u_ 

Rv'- = ~Q~  ' ^ 


sin' 


;“i  - tanh  ^)" 


(A-13) 


The  velocity  ratio  is  given  by 


j-  cosh(— 

~ sina,'  I 3 y I 

l — — “s°i"  1 - cosh(«;7r7 

n 1,1 L_ 1_J 
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(A- 14) 
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The  DC  inductance  is 


jvr  ~ A 


(sf  a*  = if 

\u  ' d2 


f+1 


J-1 


■+1 


-1 


sinaj 
i=l  “ 
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ft) 


“ a Tai 

where  d = J — r tanh  

i=i“i  La  aJ 


The  integrand  is  of  the  form 


d^dn 


(A-15) 


l yi(n)xiU) 
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2 OO  00 

= I I yi(n)xi(c)yn(n)xj(c) 
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where 


XiU) 
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cosh  — £ 
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cosh 


and 


sma^ 

yj_(n)  = 5—  cosa^ri 


The  integrand  is  sufficiently  well-behaved  to  permit  interchanging  the 
orders  of  summation  and  integration.  Thus,  the  integral  of  equation 
(A-15)  may  be  written  as 
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00 


l 

j = l 


00 
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i=l 
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yi(n)yj(n)  dn 
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/■+1 

x^(OxjU)d£ 
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(A-17 ) 
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The  integration  is  greatly  simplified  by  noting  that 


yi(n)yj<n)dn  = o,  i t j 


J-i 


i = 3 


( A-18 ) 


The  integral  is  evaluated,  and  Lvr  becomes: 


4d2 


l (z  - ~ tanh  — - tanh2  — \ 

i=l  «i6  V ai  a) 


(A-19 ) 


It  is  of  interest  to  include  the  case  of  steady  developed  laminar  flow 
between  parallel  plates,  because  this  case  represents  limiting  values 
of  both  radius  ratio  for  annular  lines  and  aspect  ratio  for  rectangular 

a 

lines.  As  r^  1 (a-*  0),  the  pairs  RVn»  Lvn  and  RVr>,  LVr  will 
approach  the  values  obtained  for  parallel  plate  flow.  Taking  a unit 
width  of  flow  between  parallel  plates  where  2b  is  the  distance  separating 
the  plates  [30 ] and  y is  measured  from  the  center  plane  between  the 
plates , \ 


■-5K  (b2-^> 


(A- 20) 


-3P/3x  , 
u = — r- b 

3 jj 


(A-21 ) 


And  the  DC  resistance  per  unit  width  is 


-9P/9x  _ 3u_ 

Q 2b3 


( A-22 ) 
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The  DC  inductance  is  easily  calculated  as 


1.2  La 


(A-23) 


Thus,  the  low  frequency  (DC)  inductance  does  not  vary  greatly  over  the 
range  of  aspect  or  radius  ratios  for  rectangular  or  annular  lines. 

With  inside  radius  r^  = 0 for  annular  lines,  Kl^  = Klc  = 4/3.  As 
r^  -*•  rc,  Ki^  -*■  1.2.  For  rectangular  lines,  has  its  maximum  at 
a = 1 where  Kj^  = 1.3785,  and  approaches  = 1.2  as  a + 0. 

A similar  approach  may  be  used  to  obtain  the  low  frequency  (DC) 
viscous  laminar  resistances  and  inductances  for  other  noncircular  lines, 
provided  the  steady  flow  velocity  profile  is  known. 
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Appendix  B 


High  Frequency  Approximations  for  Annular  Lines 


First  order  approximations  of  the  system  variables  describing  the 
performance  of  annular  fluid  transmission  lines  at  high  frequency 
ratios  (oj/wv  >>  1)  are  obtained  from  equivalent  circuit  considerations 
and  analysis  of  the  asymptotic  behavior  of  the  Bessel  function  solutions. 
In  general  at  high  frequency  ratios,  the  flow  is  approaching  inertial, 
adiabatic  behavior.  As  was  pointed  out  by  Nichols  C3],  the  velocity 
profile  is  characterized  by  a large  central  plug  wherein  inertial 
effects  are  predominant,  and  a viscous  skin  region  wherein  viscous  shear 
effects  are  predominant.  The  viscous  region  is  characterized  by  a 
skin  depth  6V  which  is  proportional  to  /v/m.  As  frequency  (w)  increases, 
<5V  decreases.  Similarly  heat  transfer  effects  may  be  separated  into  a 
relatively  large  adiabatic  core  region  bounded  by  a nonadiabatic  region, 
the  latter  being  characterized  by  a thermal  boundary  layer  thickness 
6-j>  which  is  proportional  to  /v^/us.  For  the  annular  line,  there  are  two 
boundary  layers  and  their  influence  on  the  velocity  and  temperature 
profiles  vary  with  radius  ratio.  The  behavior  is  described  in  the 
velocity  and  temperature  profiles,  equations  (11)  and  (12),  respectively. 


■W) 


has  its  greatest  influence  on 


In  general,  the  Bessel  function  J o \ 

the  solution  at  the  outer  wall  (r  = r0)  while  Ko  ^ I is  most 
influential  at  the  inner  wall.  At  high  frequencies  the  following 
inequalities  apply: 
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1 * 


J0^ro)  >:>  JC)(ri) 
Jl(r0)  >> 


K0(ri)  >>  K0(ro) 
Kj(ri)  >>  Kj(r0) 


(B-l) 


and 
energy 


J0(ro)K0(ri)  >>>  J0(ri>K0(ro)  (B-2) 

where  for  Jo»Ji  the  argument  (r)  is  presumed  to  be  (/w/v  rj3/^ 
for  KqjKj  the  argument  is  presumed  to  be  rj^^.  For  the 

equations  v is  replaced  by  v^. 

With  the  foregoing  approximations  the  parameters  appearing  in 
equations  (11)  and  (12)  may  be  greatly  simplified.  The  coefficients 
B,  D,  and  M of  equations  (13),  (14),  and  (18)  become: 


B as  B' 


*o(ri)  _ _ 1 

Jo(r,o)Ko(ri)  " " J0^ro) 


(B-3) 


Jp(ro)  _ 1 

Jo(r0)K0(ri)  “ " K0(ri) 


M as  M' 


r0J i (rQ)  . riKi (rj ) 

J0(ro)  + ^ K0(ri) 


(B-4) 

(B-5) 


From  [5,6]  the  first  order  approximations  of  the  Bessel  functions  are 
given  in  Table  B-I.  Using  the  expressions  given  in  Table  B-I  and  keeping 
only  first  order  terms,  the  high  frequency  series  impedance  becomes 


ZHI  = 3uLa  j1  + ^ ~ (r*~l)  + j(r*+l)j| 


-1 


(B-6  ) 


and  the  shunt  admittance  is 


Yhi  ■ j»ca  j 1 - (lli)  {!?  [<r*-l>  ♦ j (r*+l)J  j (B-7) 
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TABLE  B-I 


First  Order  High  Argument  Approximations  of  , Kq , Ki 


The  high  frequency  components  of  the  equivalent  circuit  are  easily 
obtained  from  equations  (B-6)  and  (B-7)  by  separating  real  and  imaginary 
parts: 


/2vui  La  1 / 

RHI  = Real(ZHI)  = - ----  = ^ ^ La  ' N 
r0(l-r») 


1-ri* 


(B-e) 


Lhi  = Imag(ZHi)/<*>  = La 


(B-9 ) 


Ghi  = Real(YHI) 


(y-1  )Ca  i 

rn(l-r*)  2 


^ {*r) 


Ca>1 


l+ri* 

l-rj* 


(B-10) 


CHI  = Imag(YHI)/w  = .Ca 


(B-ll) 
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In  the  high  frequency  regime  and  using  the  approximations  R/wL  « 1, 
G/tuC  « 1,  the  attenuation  a and  phase  lag  6 are  given  by  Weber  [11]  as 


_ LG+RC 

“HI  '1W 


= wSlC 


(B-12) 


Using  equations  (B-8)  through(B-ll)  in  the  above  equations 


*HI 


J 

1+ri*'/. 

L 

"2  Aja'“a'l0wv 

J 

1-rj*  V 

1 o ) 

2/LaC 


(B-13) 


a'-a 


®HI  = 


The  foregoing  approximations  are  accurate  to  within  1%  for  w/wv  > 100 
and  <0.9  and  are  accurate  within  10%  for  w/wv  z 10  and  r**  < 0.5. 
The  smaller  is  r^*  the  greater  is  the  accuracy  for  a given  range  of 


Appendix  C 


Derivation  of  Small  Signal  Equations  in 
Cylindrical  Coordinates 

In  this  appendix,  the  continuity,  momentum  and  energy  equations 
governing  the  small  signal  propagation  of  pressure  waves  through  a 
fluid-filled  line  are  developed  from  the  more  general  equations  of  fluid 
mechanics  and  thermodynamics.  The  equations  are  expressed  in  a cylindrical 
coordinate  system  (r,  d,  z) . The  radial,  angular  and  axial  velocities 

•f 

in  this  coordinate  system  are  ur,  ug,  and  uz,  respectively,  and  u is  the 
vector  velocity.  In  the  following  developments 

D 3 a u©  3 3 

Dt  = at  + Ur  3r  + r ae  + Uz  az 


ia  , . 

= — — ( rur) 
r ar  r 


1 3u©  5uz 
r ae~  + az” 


i a_  i_  _ai_ 

r 3r  + r2  g q2 


a2 

3z2 


Momentum  Equations 

The  Navier-Stokes  equations  of  motion  for  a compressible  fluid  are 
as  follows: 


|£ur  _ ue2j  _ 


"l£+ 

ar  1 


+ u 


Ur  2_  9u9 

r2  r2  3 e _ 


+ ±y  f^V.u) 

*1 


(C-l) 


157 


The  principal  direction  of  small  signal  wave  motion  is  the  axial 
direction;  thus,  it  is  assumed  that  the  flow  is  axially  symmetrical.  This 
further  implies  u^  = 3/36  = 0.  This  assumption  eliminates  the  momentum 
considerations  of  Equation  (C-2).  In  addition,  the  body  forces  (Fr,  Fe, 
Fz)  are  also  assumed  negligible.  The  momentum  equations  are  thus  reduced 
to  the  following: 


P 3ur  8ur 
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“Z  3z  J 

3r 
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1 3 

+ 3 y 37 
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(C-4) 

puz  3uz 

3uz") 

3£ 

1 3uz 

L TT  + Ur  3^  + 
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3z 
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r.  3r 

3z2  J 

1 3_  f 1_  3(rur)  8uz1 

3 y 3z  [r  3r  + 3z  J 


( C—  5 ) 


Th  wavelengths  of  interest  are  assumed  large  compared  to  the  cross- 
sectional  dimensions  of  the  line,  which  implies  that  the  radial  component 
of  velocity,  ur,  and  its  derivatives  may  be  neglected.  The  above 
equations  are  further  reduced  to  the  following: 


* 
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(C-6) 
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3P  + 1 u 3__  3uz 

?r  3 9r  3z 


'3u2  3u2“ 

aP  4 s2u2 

‘ 52u2 

1 3uzl 

3z  + 3 P 8z2  + W 

. 3r2 

r 3r  J 

(C-7) 


This  assumption  may  be  expressed  as  the  following  condition:  X >>  rd 
where  rd  is  the  principal  dimension  of  the  cross  section  in  the  radial 
direction.  Since  c = fX,  and  f = w/2ir  this  condition  is  equivalently- 
expressed  as 

2ttc 

* I'd  <<  (C-6) 

u u> 

where  c = speed  of  sound  in  the  line.  It  is  next  assumed  that  the 
pressure  is  uniform  across  any  line  cross  section.  This  assumption 
removes  the  dependence  of  P on  r and  eliminates  the  radial  momentum 
equation 


It  is  necessary  to  define  the  conditions  under  which  the  axial  variation 
of  u2  is  negligible,  and  to  remove  this  dependence.  The  pressure, 
velocity,  and  density  may  be  written  in  terms  of  a mean  value  and  a 
small  perturbation  about  this  value  as  follows: 

p = p + c' 

P = p + p'  (C-10) 

u = u + u' 


159 


and  the  subscript  s has  been  dropped  front  velocity,  since  only  the  axial 
velocity  remains  in  the  momentum  equation.  Substituting  these  variables 
into  equation  (C-7),  and  separating  the  steady  and  dynamic  terms  gives: 
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D'u  1 Du 
r 3r_ 


(P+p')  •jy  + (pu+pu'+Qp'+u'p' ) 


(C-ll) 
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1 Du’ 
i'  Dr 


(C-lD 


It  is  assumed  that  the  steady  part  of  the  velocity  is  constant  in  both 
the  axial  direction  and  time.  Thus  Du/Dt  = Du/Dc  = 0,  and  the  above 
equations  are  rewritten  as  follows: 
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1 Du 
r Dr 


(013) 


(p+p' ) + (P+P' )u' 


r D‘u* 

[Dr- 
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1 Du’  1 
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liquation  (01*0  can  be  further  jx-duced  by  neglecting  higher  order  terms 
involving  products  of  perturbation  quantities: 


I in' 

' Dr 


( C - 1 5 ) 


It  is  now  possible  to  solve  the  steady  flow  problem,  equation  ( C- 13). 
bv  standard  techniques.  The  combined  homogeneous  and  particular  solution 
lor  velocity  is  as  follows: 


(C-16) 


u = C,  tn  r ♦ C2  t £ £ 


where  Cj  and  C2  are  arbitrary  constants,  subject  to  the  boundary 
conditions  of  either  circular  or  annular  lines.  Tor  circular  lines,  the 
boundary  conditions  are: 


u remains  finite  at  r = 0 
u = 0 at  r = r„ 


(C-17) 


Applying  these  conditions  to  equation  (C-16)  leads  to  the  well-known 
Hagen-Poiseuille  velocity  distribution: 


3P  1 , ? 

u = - ^<r°  -r  > 

For  the  annular  line,  the  boundary’  conditions  are: 

r rjk  ~ r ’ ■ ^ ^ 


(C-16) 
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where  r*  = r/rQ 
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1 

(C-19) 


radii  which  bound  the  annular  cross  section  of  the  line.  Applying  those 
conditions  to  equation  (C-16)  gives  the  steady  developed  laminar  How 
velocity  profile  in  an  annular  region. 


(C-20) 
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Returning  to  the  dynamic  problem,  equation  (015),  it  is  necessary  to 
develop  the  conditions  for  which 


4 


3 


U 


(C-21 ) 


This  is  done  by  assuming  the  dynamic  pressure  to  be  of  the  following 
form: 


P’ 


(022) 


and  T is  the  propagation  operator  defined  in  Section  II,  which  defines 
the  lengthwise  dependence  of  the  dynamic  pressure.  Differentiating 
equation  (022)  gives 


3P' 

3z 


Apr 


e~rZ  eju,t 


(023 ) 


To  a first  approximation,  the  axial  sensitivity  of  the  perturbation 
velocity  is 


3u'  _ 1 _3PJ_ 

3z  " $ca2  3t 

Substituting  equation  (C-22)  into  (C-24),  then  differentiating 
equation  (C-24)  with  respect  to  z gives 


(C-24) 
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Comparing  equations  (C-21),  (C-23)  and  (C-25)  it  is  seen  that 


4 32u’  „ 
3 P 3z2 


3JP' 
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when 
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*4  V 


(C-26) 


For  air  at  standard  conditions  this  amounts  to  requiring  that 
w <<  3.2*107  rad/sec.  The  radian  frequencies  considered  in  this  work- 
are  well  within  this  restriction.  Using  equation  (C-26)  in  equation  (C-1S), 
and  rearranging  terms 


v 3_  r 3u’ I 3u'  _ _ 1 3P 1 1 
r 3r  r 3r  3t  p 3z 


(C-27) 


Equation  (C-27)  is  identical  with  equation  (2),  Section  II.  The 
continuity  equation  is  considered  in  the  following  development. 

Continuity  Equation . The  continuity  equation  for  a general 
compressible  fluid  expressed  in  cylindrical  coordinates  is: 


It  + r fr  (prur>  + r fe  <0U8)  + fe  (puz>  = 0 


(C- 28) 


By  previous  assumption,  ur  = ue  = 3/36  = 0;  thus 


!t  + fe(pu^  = 0 


(C-29 ) 


Expanding  equation  (C-29)  gives 


3d  3uz  3d 
3T  + P 3z~  + Uz  37  = 


(C-3C) 
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Energy  Equation.  For  a compressible  fluid  with  constant  thermal 
conductivity,  k,  the  energy  equation  is 


+ $ + k72T  - 7 . 


Dt  Dt 


at 


*r 


(C-31) 


where  h is  the  specific  enthalpy,  q is  the  internal  heat  generation 
other  than  by  viscous  dissipation,  4>  is  the  mechanical  or  viscous 
dissipation  function,  and  qr  is  the  radiation  heat  flux  vector.  For 
small  signal  propagation  the  dissipation  $ is  assumed  negligible. 
Radiation  effects  (qr)  and  internal  heat  generation  by  chemical  or 
other  like  means  (q)  are  not  considered.  With  these  assumptions, 
equation  (C-31)  reduces  to 


(C-32 ) 


Expanding  equation  (C-32)  in  cylindrical  coordinates, 
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From  the  assumptions  made  previously  on  simplifying  the  momentum 
equations,  such  as  ur  = u@  = 3/36  = 0,  the  above  equation  becomes 


(C-33) 
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1 


Since  only  the  axial  velocity  remains,  the  subscript  z is  dropped  at 
this  point.  As  before,  the  variables  are  defined  as  mean  values  and 
small  perturbation  quantities: 

P = P + p' 

P r p + p' 

T = T + T' 


These  variables  are  substituted  into  equation  (C-34)  to  obtain  the 
small  signal  energy  equation.  The  result  is 

- 3h'  3P ' , f 32T ' 1 3T'  32T'  1 (C-36) 

p 3t  “ 3t  " |_3r2  r 3r  3z2  J 

Using  the  relation  dh  = CpdT  and  dividing  by  the  quantity  pCp  gives 

3T'  1 3P ' k f32T’  1 3T'  32T'  "I  (c  ,7« 

3^“  pCp  3t  = pCp  |_3r2  r 3r  3z2  J 

The  Prandtl  number  a 2 = pCp/k,  together  with  the  relation  y/p  = v,  gives 


k 


pCp 


(C-38) 


and  vT  is  the  thermal  diffusivity.  Substituting  equation  (C-38)  into 
(C-37),  and  rearranging 
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It  is  necessary  to  determine  the  conditions  for  which 
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Assuming  the  dynamic  perturbation  of  temperature  is  of  the  following 
form 


T*  = AT  e~?z  e^Wt  (C-41) 
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then 
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(C-43  ) 


Since  the  modulus  of  T is  approximately  equal  to  w/ca  it  is  seen  by 
comparing  equations  (C-42)  and  (C-43)  that  the  condition  of 
equation  OC-40)  is  satisfied  when 

ca  ca2°2  (C-44 ) 

li)  < < - ■ — - 

V«p  V 

For  air,  o2  = .7  so  that  the  condition  expressed  in  equation  (C-44) 
is  nearly  identical  to  that  of  equation  (C-26).  The  final  form  of 
the  energy  equation  is 

VT  3 f 3T'l  3T ' BpT  gp' 

r 3r  3r  J 3t  pCp  ^t  (C-45) 

which  is  in  agreement  with  equation  (3),  Section  II. 


Temporal  Derivative  of  Density 

It  is  thought  necessary  to  show  the  derivatioi  of  equation  (22), 
since  this  equation  is  crucial  to  the  development  of  the  line  shunt 
properties.  The  density  is  assumed  to  depend  on  pressure  and 
temperature,  p = p(P,T). 

^ * (ff)T  ■*  («),, dT  lc_ 


The  isothermal  compressibility  e is  defined  as 


1 /3V\  _ 1 / 3p\ 

" V \3p/T  ' p \^3p/T 


(C-47) 


where  V is  the  specific  volume.  The  coefficient  of  thermal  expansion 


3p  is  defined  as 


_ 1 3V  . . 1 (*E\ 

pP  ‘ V 3T  " p \3TJd 


(C-48) 


For  an  ideal  gas,  6pT  = 1 and  cP  = 1. 

Equation  (C-43)  may  be  rewritten  in  terms  of  the  temporal 
derivative  as 


(C-49) 


From  the  definition  of  the  specific  enthalpy,  h: 


dh  = de  + 


(C-bO) 


where  e is  the  specific  internal  energy.  For  a perfect  gas 


CpdT  and  de  = CvdT 


(C-51) 


Using  equations  (C-51)  in  equation  (C-50)  gives 


CpdT  = CvdT  + j dp  - P 


(C-52 ) 


Rearranging  terms  and  using  equation  (C-49)  together  with  the  ideal 
gas  relations,  BpT  = eP  = 1, 


(Cp-Cv)dT  = (1-ep)  + f BpdT 
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(C-53) 


CP‘Cv  = ~pe  = Cp  i1 T. ) 


(C-54) 


The  pressure  term  in  equation  (C-49)  may  now  be  expanded  as  follows 
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Substituting  equation  (C-55)  into  equation  (C-49)  leads  directly  to 
equation  (22)  repeated  here  for  convenience: 

= ££  + «p  [ !eI  _ Hi  (22) 

at  y at  dpp  [pcp  3t  3tJ  K 1 

The  development  of  the  remaining  equations  in  Section  II  is  relatively 
straightforward. 
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